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Abs  tract 


This  study  presents  a  method  to  analyse  the  interaction 
between  flow  in  the  main  channel  and  the  flood  plain  channel 
of  a  model  river.  The  basic  concept  of  this  method  has  been 
supported  by  observation  of  the  velocity  field  and  shear 
stresses.  The  study  was  started  by  investigating  the  flow  in 
a  straight  rectangular  channel  with  smooth  and  rough 
boundaries.  The  effect  of  the  aspect  ratio  on  the 
distribution  of  the  shear  stresses  on  the  bed  and  walls  for 
both  types  of  boundaries  was  investigated.  More  precise 
values  of  the  aspect  ratio  and  related  shear  stresses  were 
obtained  for  better  definition  of  a  broad  channel,  which  was 
wide  enough  to  ensure  that  the  presence  of  walls  did  not 
interfere  with  the  main  core  of  the  flow.  A  straight  and 
symmetrical  compound  cross-section  with  rectangular  main  and 
flood  plain  channels  was  tested  and  the  results  showed  that 
momentum  is  transferred  from  the  main  channel  to  the  flood 
plain  channel  thereby  disturbing  the  velocity  field  and  bed 
shear  stress.  It  was  found  that  the  most  important 
parameters  for  slowing  down  the  flow  in  the  main  channel  and 
increasing  the  flow  in  the  flood  plain  channel  are  the 
relative  width  and  depth  of  the  main  and  flood  plain 
channels.  It  was  also  found  that  the  lateral  variation  of 
the  longitudinal  velocity  is  similar  for  various  depths  of 
flow.  A  suitable  form  of  non-dimensional  velocity 
distribution  was  introduced.  A  relationship  for  relative 
depth  and  length  scale  in  the  main  and  the  flood  plain 
channels  was  established  suggesting  that  they  are  a  measure 
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of  lateral  penetration  of  momentum  exchange  in  the  main  and 
flood  plain  channels. 

Asymmetrical  straight  channels  with  wide  main  and  flood 
plain  channels  were  tested  and  it  was  determined  that  the 
flow  pattern  did  not  differ  in  principle  from  symmetrical 
straight  channels  and  that  the  results  for  the  symmetrical 
channels  can  be  used  for  practical  purposes. 

The  analysis  was  extended  to  include  flow  in  meandering 
channels  within  a  straight  valley.  Experimental  results 
showed  that  helicoidal  motion  was  not  confined  to  the  main 
curved  channel.  The  upper  part  of  the  flow  above  the  flood 
plain  level  has  an  effective  width  approximately  equal  to 
twice  the  amplitude  of  the  meandering  channel.  In  the  main 
channel  the  flow  retained  the  secondary  flow  characteristic 
by  exhibiting  higher  velocity  along  the  radially  innermost 
(convex  wall)  wall  and  lower  velocity  along  the  outermost 
wall  (concave  wall).  The  tangential  velocity  profiles  were 
approximately  uniform  and  the  logarithmic  law  was  not  valid, 
even  as  a  rough  approximation.  The  exchange  of  momentum  was 
observed  to  take  place  between  the  channels  with  the  breadth 
of  twice  the  amplitude  of  the  main  channel  curve  and  the 
flood  plain  channel.  Comparison  of  the  flow  characteristics 
with  and  without  the  flood  plain  channel  presented  a  number 
of  interesting  points  and  showed  that  it  is  possible  to 
determine  the  zone  of  increasing  velocity,  length  of 
penetration,  and  similarity  of  lateral  velocity  profiles. 
Because  the  velocity  of  the  flow  has  a  tendency  to  persist 
in  one  direction,  it  accelerates  inside  the  curve  wave  over 
the  flood  plain  and  decelerates  in  the  main  channel.  Isovel 
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maps  showed  that  high  velocity  would  occur  on  top  of  the 
flood  plain  close  to  the  convex  region  suggesting  that  it  is 
the  excess  shear  stress  associated  with  higher  velocity 
which  erodes  the  channel.  This  mechanism  appears  to  be  the 
principal  reason  for  the  straightening-out  of  the  main 
meandering  channel  during  the  time  of  high  flood  in  a  river. 

Qualitative  tests  were  also  conducted  in  the  straight 
and  meandering  channels  using  dye  to  observe  the  flow 
direction.  Eddy  viscosity  computations  were  made  in  a  few 
cases  . 


v  i 


ACKNOWLEDGMENTS 


The  author  wishes  to  express  his  sincere  thanks  to  Dr.  N. 
Rajaratnam  for  his  guidance,  advice  and  encouragement 
thoughout  the  course  of  this  work.  Thanks  are  due  to  Dr. 

J.P.  Verschuren  for  his  interest  in  this  work,  to  Dr.  G. 
Parker  for  many  intereting  discussions  on  this  problem,  Dr. 

R.  Torgerson  for  his  advice  and  guidance  to  adopt  the 
Surface  II  Program  in  preparation  of  the  isovel  maps  and  to 

S.  Lovell,  D.  McGowan,  R.  Gitzel  and  A.  Muir  for  their 
cooperation  and  help  in  the  construction  of  the  experimental 
arrangement  and  in  the  instrumentation. 

The  author  is  grateful  to  the  University  of  Alberta  for 
providing  a  dissertation  fellowship  and  to  the  National 
Pesearch  Council  of  Canada  for  partial  financial  support  of 
this  study  through  a  grant  to  Dr.  N.  Rajaratnam. 


vii 


. 


, 


Table  of  Contents 


Chapter  ,  Page 

1.  Introduction . 1 

1  .  1  General . 1 

1.1.1  Introduction  to  Problem . 1 

1.2  Review  of  Flow  in  Channels  with  Flood  Plain . 2 

1.3  Effect  of  Curvature  of  Fain  Channel . 3 


1.3.1  Review  of  the  Literature  on  Curved  Channels.  ..4 

1.3.2  Comments  on  Assumptions  in  Flow  Around 


Bends . 6 

1.4  Scope  of  Present  Study . 10 

1.4.1  Presentation  of  the  Program . 11 

2.  Experimental  Arrangement  and  Measurement  Techniques  ....  1 5 

2.1  Introdution . 15 

2.2  Experimental  Flume . 15 

2.2.1  Boundary  Roughnes . 16 

2.2.2  Temperature  Measurements . 16 

2.3  Experimental  Channel . 16 

2.3.1  Symmetrical  Channel . 17 

2.3.2  Asymmetrical  Channel . 18 

2. 3. 2.1  Asymmetrical  Channel  Type  A . 18 

2. 3. 2. 2  Asymmetrical  Channel  Type  B . 20 

2. 3. 2. 3  Asymmetrical  Channel  Type  C . 20 

2. 3. 2. 4  Asymmetrical  Channel  Type  D . 20 

2.3.3  Meandering  Channel . 21 

2.4  Velocity  Measurements . 24 

2.5  Velocity  Vector  Measurement . 25 


viii 


. 

.  .  . 

•  I 


2.6  Boundary  Shear  Stress  Measurements . 29 

2.6.1  Yaw  Probe  as  a  Preston  Tube . 30 

2.7  Depth  and  Slope  of  Water  Surface . 32 

2.8  Special  Instrumentation . 33 

2.8.1  Coordinate  Positions  and  Input  Panel . 33 

2.8.2  Data  Acquisition  System . 35 

2.8.3  Test  Method . 37 

2.9  Flow  Visualisation . 38 

3.  Interaction  Between  Straight  Main  Channel  And  Flood 

Plain  Flows  in  a  Symmetrical  Channel . 39 

3.1  Concept  of  Momentum  Exchange  and  Interaction 

Region . 39 

3.2  Experiments  And  Experimental  Results . 40 

3.3  Analysis  of  Experimental  Results . 46 

3.4  Effect  of  Roughness . 62 

3.5  Conclusions . 62 

4.  Interaction  Between  Straight  Main  Channel  and 
Straight  Flood  Plain  Flows-Further  Studies-Asymmetri- 

c  al  Ch  annel . 65 

4.1  Experiments . 65 

4.2  Experimental  Results . 66 

4.3  Analysis  of  Results . 74 

4.4  Kinematic  Eddy  Viscosity . 93 

4.5  Conclusions . 98 

5.  Meandering  Main  Channel  with  Straight  Flood  Pla  ins  .  .  .  .  1 00 

5.1  Flow  in  Curved  Channels . 100 

5.2  Theoretical  Considerations  of  Flow  in  a  Bend . 102 

5.3  Experiment  and  Experimental  Results  on  Flow  in 

Curved  Channels . 105 

ix 


5.4  Experiments  on  Meandering  Channels  with  Flood 

Plains . . 

5.5  Flow  Visualization  Results . 117 

5.6  Data  Handling  and  Presentation  of  Pesults . 118 

5.6.1  Velocity  Distribution . 120 

5.6.2  Cross-Sectional  Isovel . 179 

5.6.3  Isovel  Map . 184 

5.7  Discussion  on  Interaction  Results . 192 

5.8  Conclusions  and  Comments  . 194 

5.9  Further  Direction . 195 

6.  Summary  of  Conclusions . 197 

6.1  Concept  of  Interaction  between  Main  Channel  and 

Flood  Plain  Flows . 197 

6.2  Hide  Main  Channel  and  Hide  Flood  Plain  Channel. ...197 

6.3  Meandering  Channel  in  a  Straight  Valley . 199 

NOTATIONS . 201 

References . 205 

AP PND I X  A . 212 

A.-1  Turbulent  Flow  in  Rectangular  Channels . 212 

A-2  Experimental  Arrangements  and  Experiments . 213 

A-3  Velocity  Distribution . 214 

A -4  Distribution  of  Boundary  Shear  Stresses . 215 

A- 5  Average  Velocity  Prediction . 217 

APPENDIX  B . 220 

B-l  Channel  with  Flood  Plains  at  Two  Levels . 220 

B-2  Analysis  of  Experimental  Results . 225 

B-3  Conclusions . 226 


x 


' 


Tables 

T  ab le  Page 

Table  1  Significant  Details  of  Experiments 

in  Chapter  III .  41 

Table  2  Significant  Details  of  Experiments 

in  Chapter  IV . 67 


Table  of  Figures 


Figure 
Figure  1 
Figure  1 
Figure  2 
Figure  2 

Figure  2 
Figure  2 

Figure  2 

Figure  2 

Figur  e  2 

Figure  2 
Figure  2 
Figure  3 
Figure  3 

Figure  3 

Figure  3 
Figure  3 

Figure  3 
Figure  3 


1 

2 

1 

1A, 

2 

3 

4 

5 

6 

7 

8 
1 

2a, 

3a  , 

4  a 

5 

6 
7 


Page 


Variation  of  Bed  Slope . 

Flow  Around  Bluff  Body . 

Geometry  of  Straight  Channels . 

B  &  C  Geometry  of  Asymmetrical  Channel 

Types  A,  B,  C  and  D . 

Geometry  of  Curved  Channel . 

Location  of  Cross-Sections  in  Curved 

Channel . 

Yaw  Probe  Calibration  curve  (Variation 

of  k^,  k£  &  k^  with  9) . 

Yaw  Probe  Calibration  curve  (Variation 

of  k  with  9) . 

Yaw  Probe  Calibration  curve  (Variation 

of  k  with  9) . 

o 

Photograph  of  Data  Acquisition  System.... 

Flow  Chart  of  Equipments . 

Defination  Sketches . 

b  &  c  Typical  Velocity  Profiles  (in  y  &  y' 

direction) . 

b  &  c  Typical  Velocity  Profiles  (in  z 

direction) . 

&  b  Bed  Shear  Stress  Profiles . 

Velocity  Profile  Similarity:  Upper  Part 

of  Main  Channel . 

Velocity  Scale:  Upper  Part  Main  Channel.. 
Length  Scale  :  Upper  part  of  Main  Channel 


7 
7 
1  9 

1  9 
22 

22 

28 

28 

28 

34 

36 

43 

44 

45 
47 

49 
5  1 
52 


x  i  i 


. 


t 


54 


Figure  3.8  Velocity  Profile  Similarity:  Flood  Plain. 

Figure  3.9  Velocity  and  Length  scales:  Flood  Plain..  55 

Figure  3.10  Velocity  Profile  Similarity:  Lower  Part 

of  Main  Channel .  56 

Figure  3.11  Bed  Shear  for  Flood  Plain  and  Defect  in 

Main  Channel .  58 

Figure  3.12  Centerline  Bed  Shear  for  Rectangular 

Channels .  60 

Figure  4.1  Definition  Sketches .  68 

Figure  4.2  Velocity  Profile:  Expt-3 .  69 

Figure  4.3  Lateral  Velocity  Profiles:  Expt-1,  3  &  4.  71 

Figure  4.4  Bed  Shear  Stress:  Expts  1  to  8 .  75 

Figure  4.5  Velocity  Profile  Similarity:  Upper  Part  of 

Main  Channel  (Expts  1  to  7) .  78 

Figure  4.6  Calculated  Velocity  vs.  Measured  Velocity  81 

Figure  4.7  Variation  of  Normalized  Depth  versus 

Normalized  Velocity .  82 

Figure  4.8  Length  Scales:  b^,  b^,  b  &  bT  .  83 

Figure  4.9  Velocity  Profile  Similarity:  Flood  Plain 

(Expts  1  to  6) .  8  5 

Figure  4.10  Bed  Shear  Stress  Similarity:  Flood  Plain 

Channel  (Expts  1  to  8) .  87 

Figure  4.11  Variation  of  Bed  Shear  Stress  with  D/d  &  .89 

Figure  4.12  Velocity  Distribution  (u(y')  & 

u  (y  )  )  :  l  og  Law .  90 

Figure  4.13  Location  of  Dip  in  Velocity  Profiles:  Main 

Channel .  94 

Figure  4.14  Eddy  Viscosity . .  96 

Figure  4.15  Cross-Sectional  Isovel  Map:  Fxpt-4 .  96 


xi  ii 


■  . 

I 


107 


Figure  5.1  Distribution  of  Tangential  Velocity  in 

Curved  Channel,  Cross-Sections  1  to  15... 
Figure  5.2  Distribution  of  Fa  dial  Velocity  in 

Curved  Channel,  Cross-Sections  1  to  15...  108 
Figure  5.3  Photograph:  Helicoidal  Motion  of 

the  Flow  in  the  Meandering  Channel .  110 

Figure  5.4  Cross-Sectional  Isovel  Map  in  Curved 

Channel  (Cross-Sections  1  to  15) .  Ill 

Figure  5.5  Lateral  Distribution  of  u  . 112 

Figure  5.6  Bed  Shear  Stress  Distribution .  112 

Figure  5.7  Variation  of  Velocity  and  Depth  Along 

Centerline . 113 

Figure  5.8  Non-dimensional  Centerline  Velocity 

Profiles  in  Curved  Channel,  Cross-Sections 

1  to  15 .  113 

Figure  5.9  Non  Dimensional  Lateral  Distribution  of 

Tangential  Velocity  in  Curved  Channel, 

Cross-Sections  1  to  15 .  113 

Figure  5.10  Photograph:  Efective  Zone  of  Momentum 

Exchange  in  the  Meandering  Channel 

with  Flood  Plain  Channels  .  119 

Figure  5.11(1  -  15)  Angle  of  Velocity  Vector  Profiles: 

Cross-Sections  No .  1  to  15 .  121 

Figure  5.12(1  -  15)  Padial  Velocity  Profiles: 

Cross-Sections  No.l  to  15 .  133 

Figure  5.13(1  -  15  )  Velocity  Vector  Profiles: 

Cross-Sections  No.l  to  15 .  145 

Figure  5.14(1  -  15)  Tangential  Velocity  Profiles: 

Cross-Sections  No.l  to  15 .  157 


xiv 


Figure  5.15  Lateral  Tangential  Velocity 

Distribution:  Cross-Section  No.  0  to  15 

CP-UNl) .  172 

Figure  5.16  Lateral  Tangential  Velocity 

Distribution:  Cross-Section  No.  1  to  15 

C RUN 2  ) .  178 

Figure  5.17  Cross-Sectional  Isovel  Nap:  RUN1 .  180 

Figure  5.18  Tangential  Velocity  Isovel  Map  at 

Different  Level  of  y '  :  RUNl .  185 

Figure  5.19  Angle  of  Velocity  Vector  Contour  Map  at 

y ' =0 . 0  0  5  Feet:  RUNl .  187 

Figure  5.20  Velocity  Vector  Isovel  Map  at 

Different  Level  of  y'  :  RUNl .  188 

Figure  5.21  u^,  V  &  0  Contour  Map  at  y ' =  0 . 0 5  Feet: 

RUN  2 .  189 

Figure  5.22  Radial  Velocity  Isovel  Map:  RUNl  £•  RUN 2..  190 

Figure  A-l  Centerline  Bed  Shear  Stress  for 

Rectangular  Channels:  Smooth  &  Rough 

Boundary .  216 

Figure  A -2  Lateral  Distribution  of  Bed  Shear  Stress 

for  Pectangular  Channels:  Smooth 

Boundary .  216 

Figure  A-3  Lateral  Distribution  of  Bed  Shear  Stress 

for  Rectangular  Channels:  Rough 

Boundary .  216 

Figure  B-l  Velocity  Profiles:  Main  Channel,  First 

&  Second  Flood  Plain  Level  Channel .  221 

Figure  B-2  Lateral  Velocity  &  Bed  Shear  Stress 

Distribution .  224 


xv 


. 


♦ 


1 .  Introduc t ion 


1  .  1  General 

If  high  flow  in  a  river  rises  above  the  natural  or 
artificial  banks  in  any  reach,  water  from  the  overtopped 
banks  spreads  out  over  the  flood  plain.  Since  the 
characteristics  of  the  flood  plain  make  it  a  desirable 
location  for  a  variety  of  activities,  it  is  important  to 
develop  a  method  to  determine  the  rating  curve  for  high 
discharges . 

1.1.1  Introduction  t  o  Problem 

Measures  to  control  flood  flows  are  only  one  set  of 
adjustments  that  can  be  made  to  combat  the  flood  hazard.  The 
appraisal  of  such  engineering  measures,  if  they  are  to  be 
viewed  in  relation  to  their  full  consequences,  requires  the 
analysis  of  other  possible  adjustments.  Experience  shows 
that  the  construction  and  operation  of  engineering  works  to 
control  flood  flows  may  reduce  or  curb  flows  without 
necessarily  reducing  the  threat  of  erosion  to  the  river  bed 
and  banks  or  the  flood  plain  area.  Little  is  known  of  the 
shear  stress  distribution  on  a  river  bed  and  flood  plain 
bed.  Rise  of  flow  over  the  flood  plain  area  together  with 
the  increase  of  shear  stress  is  capable  of  eroding  the  edge 
of  the  flood  plain  and  the  main  channel.  The  redistribution 
of  velocity  and  shear  stress  will  affect  the  amount  of 
sediment  transported  in  both  channels  as  well  as  the 
stability  of  the  main  channel  bank  and  the  flood  plain  beds. 
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1  »  2  ?.  e  v lew  o f  F 1  ov  i  n  Channels  v i  th  Flood  Plain 

When  the  depth  of  flow  D  in  the  main  channel  exceeds 
the  depth  of  the  main  channel  h,  the  flood  plain  carries  a 
part  of  the  total  discharge.  In  the  preparation  of  rating 
curves  for  rivers  when  D  is  much  larger  than  d  (where  d  is 
the  depth  of  flow  in  the  flood  plain)  the  general  practice 
(Chow,  1959)  is  to  divide  the  compound  flow  section  into  a 
number  of  homogeneous  sub-areas  by  the  introduction  of 
shear-free  vertical  boundaries  which  are  the  extension  of 
the  banks  of  the  rain  channel.  The  flow  is  then  computed  for 
each  subsection  and  the  s ur  of  the  flow  in  all  of  the 
subsections  is  the  total  discharge  0.  When  D  is  only 
slightly  larger  than  d,  the  division  is  done  by  introducing 
a  horizontal  plane  as  an  extension  of  the  bed  of  the  flood 
plain  (Smith,  197F,  and  Toebes  &  Sooky,  1967).  The  first 
method  of  sub-division  is  the  one  most  frecuentlv  used  since 
the  majority  of  the  practical  cases  have  a  large  value  of 
D/d.  It  has  been  generally  recognized  (Goncharov,  1964, 
Sellin,  1964,  and  Zhe le zny akov ,  1971)  that  appreciable 

interaction  and  turbulent  mixing  occurs  between  the  faster 
main  channel  flow  and  the  slower  flood  plain  flow.  Because 
of  the  transport  of  the  longitudinal  momentum  to  the  flood 
plain,  the  vertical  plane  between  the  main  ch  annel  and  flood 
plain  experiences  significant  (turbulent)  shear  stress  in 
the  longitudinal  (or  flow)  direction.  Vright  &  Carstens 
(1970)  demonstrated  the  existence  of  this  shear  stress  in  an 
air-duct  of  a  c omp ound  cross-section.  A  similar  observation 
could  be  made  from  the  study  of  byers  &  Elsawv  (1975)  on  a 
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rectangular  main  channel  with  an  asymmetrical  flood  plain. 

Sellin  (1964)  found  that  when  the  depth  of  flow  D  just 
exceeds  the  depth  of  the  main  channel  h,  the  rate  of 
increase  of  the  discharge  Q  with  the  total  depth  of  flow  D 
becomes  very  small.  For  cases  with  a  large  value  of  D/d  the 
discharge  Q  of  the  compound  cross-section  could  be 
considerably  less  than  the  sum  of  the  individual  discharges 
of  the  homogeneous  sub-areas,  obtained  by  the  introduction 
of  the  shear-free  vertical  separating  walls.  A  few 
researchers  (Sellin,  1964,  Zheleznyakov ,  1971,  and  Myers  & 

Elsawy,  1975)  found  that  in  a  channel  with  an  asymmetric 
flood  plain,  the  bed  shear  stress  in  the  flood  plain  would 
be  much  larger  than  its  undisturbed  value,  which  would  exist 
if  there  were  no  interaction  with  the  main  channel  flow. 
Similarly,  the  main  channel  bed  would  experience  a  reduction 
of  the  bed  shear  stress  from  its  undisturbed  value.  Some 
further  observations  on  the  interaction  in  straight  channel 
have  been  made  by  Sellin  (1964),  Zheleznyakov  (1971)  and 
Ghosh  &  Jena  (1971). 


1 . 3  Effect  o  f  Curvature  o  f  Main  Channel 

It  is  well  known  that  natural  channels  are  neither 
straight  and  symmetrical  nor  are  the  boundaries 
hydraulically  smooth.  L.B.  Leopold  and  W.B.  Langbein  in  a 
United  States  Geological  Survey  report,  "River  Meanders" 


(  1  9  6  6)  wr o  t  e : 
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"Is  there  such  a  thing  as  a  straight  river?  Almost 
anyone  can  think  of  a  river  that  is  more  or  less 
straight  for  a  certain  distance,  but  it  is  unlikely 
that  the  straight  portion  is  either  very  straight  or 
very  long.  In  fact,  it  is  almost  certain  that  the 
distance  any  river  is  straight  does  not  exceed  10 
times  its  width  at  that  point." 


1.3.1  Peview  o  f  the  Literature  o  n  Curved  Channels 

Thomson  in  1876  observed  the  spiral  motion  in  a  river 
bend.  Boussinesq  theoretically  solved  the  equations  of 
motion  for  the  case  of  laminar  flow  in  a  curved  pipe.  Many 
investigations  have  since  been  conducted  on  flow  in  bends. 
P.C.  Yen  (1965)  gives  a  brief  summary  of  the  important 
experiments  on  flow  in  open-channel  bends. 

Fargue  summarized  his  observation  of  the  Garonne  river 
since  1849  and  proposed  a  number  of  empirical  laws  dealing 
with  a  meandering  river  and  a  mobile  bed.  Leliavsky  verified 
Fargue's  laws  from  field  observations.  In  1934,  Elue, 
Herbert,  and  Lancefield  published  results  of  measurements  in 
the  Iowa  Fiver.  In  1935,  Eakin  published  a  study  on  the 
curvature  of  the  Mississippi  Fiver. 

Among  the  analytical  and  experimental  investigators, 
Boss  assumed  free-vortex  velocity  distribution  in  the  radial 
direction.  Yarnell  and  Woodward  observed  spiral  motion. 
Mockmore  performed  a  theoretical  study  on  flow  in 
open-channel  bends  by  using  assumed  distributions  of 
velocity  components.  His  study  showed  the  streamlines  to  be 
helicoidal. 

Shukry  (1949)  presented  the  results  of  his  experiments, 
which  covered  a  wide  range  of  central  angle,  Froude  number. 
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and  aspect  ratio. 

Rozovskii  C  1  9  5  7  )  reviewed  the  work  of  previous  Russian 
investigators  such  as  An any an  and  Milovich.  With  the  aid  of 
o r d e r-o f -ma gn i t ud e  considerations  and  assumptions  about  eddy 
viscosity,  vertical  distribution  of  longitudinal  velocity 
components,  and  zero  net  lateral  discharge,  Rozovskii  was 
able  to  derive  an  approximate  solution  of  the  radial 
velocity  component  from  Reynolds  equations  of  motion.  Using 
a  number  of  other  assumptions  he  attempted  to  solve  f  o  r  the 
radial  distribution  of  the  longitudinal  velocity  component 
and  for  the  growth  and  decay  of  the  spiral  motion.  The 
classical  work  in  this  area  was  done  by  Rozkovskii.  Wore 
contemporary  investigators  have  tended  to  accept  Pozovskii's 
a  s  sump  t i o n s . 

A  widely  accepted  analytical  and  experimental 
investigation  is  that  conducted  by  R.C.  Yen  in  1965.  Yen 
proposed  that,  for  fully  developed  flow  in  a  bend,  the  loss 
of  energy  is  the  same  for  equal  increments  of  the  turning 
angle.  Based  on  a  logarithmic  distribution  of  longitudinal 
velocity  components  along  the  vertical,  and  assuming  that 
the  eddy  viscosity  is  a  function  of  the  depth  and  the  radial 
distance,  he  derived  an  expression  for  the  vertical 
distribution  of  the  radial  velocity  component.  However,  by 
assuming  that  the  coefficient  of  roughness  in  a  curved 
channel  is  the  same  as  that  in  a  straight  channel,  he  was 
able  to  derive  an  expression  for  the  distribution  of  the 
average  longitudinal  velocity  in  the  radial  direction.  He 
found  that  the  average  longitudinal  velocity  is 
approximately  inversely  proportional  to  the  square  root  of 
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the  radial  distance  from  the  center  of  the  meander  curve. 

In  1974  Fngelund  assumed  a  constant  tangential  slope 
for  fully  developed  flow  in  the  curved  channel  and  used 
similarity  assumptions  for  radial  and  longitudinal  velocity 
components.  Further,  he  assumed  that  the  longitudinal 
velocity  followed  the  logarithmic  distribution  near  the  bed 
and  parabolic  distribution  in  the  upper  part  of  the  flow. 
Eased  on  these  considerations,  Fngelund  derived  a  fairly 
simple  equations  for  the  distribution  of  radial  velocity 
over  the  depth  and  of  longitudinal  velocity  over  radial 
distance . 

Ikeda's  (1975)  work  closely  followed  that  of  Fngelund. 
However,  Ikeda  proposed  an  arbitrary  function  for  the 
distribution  of  tangential  velocity  over  the  width  of  the 
channel  . 

1.3.2  Comment  s  o n  A s  sump  t ions  i n  Flow  Around  Bends 

In  order  to  simplify  the  Reynolds  equations  of  motion, 
certain  general  assumptions  have  been  made  by  all 
investigators.  An  assumption  which  is  frequently  employed  is 
the  idea  of  fully  developed  flow  around  a  bend.  The  fully 
developed  flow  in  a  pipe  is  reported  after  300  degrees  turn 
or  more  (Ellis  &  Joubert,  1974).  In  considering  the  flow  in 
an  open  channel,  it  is  impossible  to  reach  fully  developed 
flow  because  of  the  presence  of  the  free  water  surface.  If 
the  bed  slope  of  the  valley  is  assumed  constant  throughout 
the  entire  reach  as  one  can  conclude  from  Figure  1.1,  the 
bed  slope  of  the  main  channel  (in  the  straight  portion)  is 
the  same  as  the  valley  slope.  But  as  one  moves  along  the 
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Magnitude,  of  Bed  Slope 


Figure  1.1  Variation  of  Bed  Slope 


FI  ow 


Figure  1.2  Flow  Around  Bluff  Body 
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meandering  path,  the  bed  slope  of  the  main  channel  changes 
accordingly.  When  the  channel  is  running  cross- wise,  the  bed 
slope  is  zero.  If  one  assumes  a  flat  bed,  a  water  surface 
slope  is  required  in  order  that  the  flow  can  take  place,  and 
since  this  means  a  drop  in  elevation  the  depth  of  flow  will 
not  be  the  same  at  all  cross-sections.  Thus,  one  cannot 
assume  there  is  uniform  flow,  with  a  constant  slope  and 
depth  of  flow,  even  at  the  centerplane  of  the  meandering 
channel.  Choudhary  et  al.  (1977)  observed  that  in  a  180 
degree  bend,  of  a  narrow  rectangular  channel  (aspect  ratio 
of  5),  and  of  a  wide  rectangular  channel  (aspect  ratio  of 
20),  the  depth  of  flow  changes  with  the  angle  of 
cross-section  0.  Sieber  &  Gotz  (1975)  investigated  the 
deviation  of  the  tangential  velocity  from  the  mean  velocity 
with  0.  B.C.  Yen  and  C.L.  Yen  (1971)  studied  the  variation 
of  the  water  surface  slope  with  0. 

Another  assumption  is  R=S  r  in  which  Sq  represents 
the  water  surface  slope  along  the  channel  centerline,  R  is 
radius  of  thee  centerline  and  S  is  the  water  surface  slope 
at  a  radial  distance  of  r  from  the  center  of  the  bend.  B.C. 
Yen  (1965)  arrives  at  this  conclusion  by  setting  the  loss  of 
energy  equal  for  every  segment  of  egual  included  angle. 
Further, it  is  well  known  that  flow  will  not  behave  the  same 
before  and  after  it  passes  a  bluff  body  (see  Figure  1.2). 

Even  with  these  assumptions,  the  number  of  unknowns  in 
the  simplified  equations  of  motion  are  still  more  than  the 
number  of  equations  of  motion  and  continuity.  To  circumvent 
this  difficulty,  many  investigators  have  had  to  assume  a 
vertical  distribution  of  the  tangential  velocity.  With  few 
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exceptions,  investigators  have  assumed  that  the  logarithmic 
law  is  adequate  to  describe  the  vertical  distribution  of 
this  tangential  velocity  component.  Rozovskii,  unlike  the 
others,  considered  the  parabolic  law,  the  elliptic  law  and 
the  logarithmic  law,  but  preferred  to  use  the  logarithmic 
law  with  the  Von  Karman  constant  equal  to  0.5.  If  the 
momentum  exchange  can  occur  in  the  vertical  direction  from 
the  upper  layer  to  the  lower  layer,  there  is  no  reason  to 
believe  that  a  momentum  exchange  can  not  occur  in  horizontal 
direction  from  the  inner  region  (generally  a  region  of 
higher  velocity)  toward  the  outer  region  (the  region  of 
lower  velocity).  It  is  known  that  v ' ,  the  vertical 
fluctuation  velocity,  which  is  only  a  small  percentage  of 
the  tangential  velocity  u,  is  responsible  for  the  vertical 

v7 

mixing.  Although  the  magnitude  of  the  radial  velocity  v  , 
depends  upon  many  factors  such  as  the  radius  of  curvature, 
the  aspect  ratio,  its  order  of  magnitude  is  comparable  to 
that  of  v'.  Therefore,  the  tangential  velocity  must  be 
redistributed  in  the  lateral  direction.  As  far  as  this 
author  has  been  able  to  determine,  no  one  has  accounted  for 

the  effect  of  v  on  the  redistribution  of  urt,  over  either 

r  9 

the  depth  or  the  width  of  the  flow. 

Similarity  assumptions  for  the  distribution  of  u^  and 
v  appear  to  be  unfounded.  In  sections  beyond  the  crest  of 
the  bend  (i.e.  in  a  180  degree  bend,  sections  beyond  90 
degrees),  the  trend  of  the  predicted  lateral  distribution  of 
uQ  is  contrary  to  that  existing  in  the  flow  (R.  C.  Yen, 

1965,  Choudhary,  1977,  and  Siebert  &  Gotz,  1975). 

An  adequate  model  should  not  only  be  able  to  predict 
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the  distribution  of  radial  velocity,  but  should  also  be  able 
to  account  for  the  variation  with  0  as  well  as  the  effect  of 
the  secondary  flow  on  the  redistribution  of  the  velocity 
field.  To  date,  investigators  are  uncertain  as  to  whether 
the  kinematic  eddy  viscosities  (£  e  e  )  in  a  straight 

X  *  y  ’  2  0 

channel  are  the  same  as  those  (£  £  ,  £  )  in  a  meandering 

w  r  z 

channel.  It  is  generally  believed  that  secondary  currents  in 
a  meandering  channel  are  stronger  than  those  in  a  similar 
straight  channel.  Therefore,  there  is  no  basis  for  assuming 
that  the  eddy  viscosities  are  even  approximately  the  same  in 
straight  and  meandering  channels. 


1  . 4  Scope  o  f  Present  Study 

An  understanding  of  the  velocity  field  and  shear  stress 
distribution  in  the  main  and  flood  plain  channels  will  be  of 
great  practical  use,  primarily  for  the  protection  of  flood 
control  structures.  Erosion  of  the  bed,  river  bank  and 
surrounding  flood  area,  conveyance  of  flow  through  the 
compound  cross-section  channel,  magnitude  of  shear  stress 
and  eddy  viscosity,  both  horizontal  and  vertical  are 
required  to  predict  the  flood  stage,  erosion  and  deposition, 
and  the  amount  of  sediment  transported  in  the  system. 

Another  factor  is  the  momentum  exchange  between  the  main 
channel  and  flood  plain  channel.  Even  though  a  certain 
amount  of  work  has  already  been  done  on  some  aspects  of  the 
problem,  many  questions  remain  to  be  answered. 

The  present  investigation  was  undertaken  in  an  effort 
to  understand  such  characteristics  of  flow  as  the 
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distibution  of  velocity,  the  bed  shear  stress,  the  effective 
zone  of  momentum  exchange  and  the  variation  of  eddy 
viscosity.  These,  if  properly  understood,  would  enable  the 
hydrologist  to  design  flood  control  structures  and  predict 
the  flood  hazard  zone  with  some  degree  of  confidence. 

The  experiments  began  with  the  simplest  possible  case: 
a  straight  rectangular  main  channel  with  symmetrical  flood 
plain.  In  an  attempt  to  study  a  main  channel  of  large  width, 
a n  a symme trical  channel  with  larger  main  channel  width  was 
tested.  The  study  was  then  extended  to  determine  whether  a 
one  stage  flood  plain  could  be  superimposed  on  another  to 
provide  a  multi-stage  flood  plain. 

Natural  rivers,  of  course,  differ  significantly  from 
the  ideal  case  considered  above.  In  nature,  rivers  are  not 
straight,  and  the  main  channel  meanders  along  its  valley  in 
a  somewhat  irregular  manner.  In  order  to  assess  the  effect 
of  the  meandering  main  channel,  the  present  investigation 
includes  the  case  of  a  curved  main  channel  in  a  straight 
valley.  Some  dye  injection  experiments  were  performed  to 
explore  the  aspect  of  secondary  flow  in  the  main  channel  and 
the  encroachment  of  secondary  flow  into  the  flood  plain 
region.  In  the  course  of  this  study,  some  supplementary 
experiments  were  performed  to  study  the  variation  of  the  bed 
shear  stress  with  the  aspect  ratio  for  simple  rectangular 
channels  with  smooth  and  rough  beds. 

1.4.1  Pr es  ent a t ion  o  f  the  Program 

The  experimental  arrangements,  the  measuring 
techniques,  the  structure  of  the  base  flume  and  straight  and 
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meandering  channels  are  discussed  in  Chapter  II.  This 
chapter  also  contains  a  description  of  the  probes,  the 
calibration  techniques  and  the  method  of  using  the 
calibration  charts.  The  limitations  of  the  measurements  are 
developed  fully  in  Chapter  II  as  well.  An  x-y-z  electric 
coordinate  system,  in  conjunction  with  a  data  acquisition 
system  were  used  to  handle  the  mass  of  data  obtained  from 
experiments.  The  data  could  be  processed  by  computer  with 
great  ease.  The  system  is  described  without  going  into 
details  of  the  electronic  design  or  programming  of  the 
mi n i-c omp u t e r  (PDP8)  Gr  the  main  computer.  A  dye  injector 
was  constucted  for  flow  visualisation  tests. 

In  Chapter  III,  investigations  on  the  interaction 
between  a  straight  main  channel  and  symmetrical  flood  plain 
flows  with  smooth  and  rough  beds  are  presented.  This 
preliminary  study  clarified  the  concept  of  momentum  exchange 
between  the  main  channel  and  the  flood  plain  flows.  The 
interaction  region  and  reduction  of  shear  stress  in  the  main 
channel  and  increase  of  the  bed  shear  stress  at  the  edge  of 
the  flood  plain  are  also  studied. 

Investigations  of  the  effect  of  a  wide  main  channel  and 
large  width  flood  plain  channel,  are  presented  in  Chapter 
IV.  Experimental  findings  on  the  asymmetrical  channel 
confirmed  the  basic  concepts  established  in  Chapter  III. 

Once  again,  it  was  possible  to  show  that  the  lateral 
velocity  profiles  are  similar.  The  length  scales  in  the 
interaction  region  were  correlated  with  D/d.  Redistribution 
of  shear  stress  in  the  main  and  flood  plain  channels  was 
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In  Chapter  V,  a  preliminary  investigation  of  the  flow 
in  a  curved  channel  without  the  flood  plain  is  presented. 

The  experimental  results  of  this  run  indicated  the 
inadequacy  of  the  existing  solutions  to  predict  the  velocity 
profiles  and  the  bed  shear  stress  in  curved  open  channels. 
Therefore,  the  results  of  this  study  on  curved  channels  are 
presented  in  an  empirical  form.  This  chapter  also  presents 
the  results  of  two  experiments  with  a  curved  main  channel  in 
a  straight  valley. 

A  number  of  point  velocities  were  measured  with  the  use 
of  the  data  acquisition  system  in  order  that  the  data  could 
be  fed  to  the  computer.  After  plotting  the  velocity  profiles 
with  the  computer  plotter,  for  better  visualization  of 
secondary  flow  and  acceleration  and  deceleration  of  flow, 
isovel  maps  in  the  y-z  and  x-z  planes  were  prepared.  These 
velocity  profiles  and  isovel  maps  showed  the  interaction 
between  the  main  meandering  channel  and  the  flood  plain 
ch  annel . 

Based  on  the  present  studies,  general  conclusions  are 
presented  in  Chapter  VI.  A  summary  of  the  idea  of 
interaction  in  a  wide  main  channel  and  wide  flood  plain,  and 
also,  a  general  qualitative  study  of  flow  in  the  meandering 
channel  in  a  straight  river  valley  are  discussed.  Some 
recommendations  are  then  made  for  further  studies  in  this 
f  ield  . 

During  the  course  of  this  investigation,  some 
supplementary  experiments  were  conducted  in  rectangular 
channels  with  smooth  and  rough  beds.  The  results  of  the 
correlation  between  maximum  bed  shear  stress  in  a 
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rectangular  channel  with  aspect  ratio  and  distribution  of 
bed  shear  stresses  are  presented  in  Appendix  A.  Appendix  B 
presents  the  results  of  one  experiment  in  a  straight  channel 
with  two  asymmetrical  flood  plains  at  two  different  levels. 

The  present  experimental  investigations  were  carried 
out  at  the  T.  Blench  hydraulics  laboratory  of  the  University 
of  A.lberta,  Edmonton,  Alberta,  Canada,  during  the  years  1977 


and  1978 
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2.  Experimental  Arrangement  and  Measurement  Techniques 


2 . 1  Introdution 

This  chapter  describes  the  construction  of  the  flume, 
the  experimental  arrangements  used,  the  limitation  of  the 
measurements  and  physical  limitation  of  the  size  of  the 
channels.  The  first  set  of  experiments  were  conducted  in  a 
main  channel  with  symmetrical  flood  plain.  The  second  set  of 
experiments  were  conducted  in  a  straight  but  asymmetrical 
channel  to  simulate  a  very  wide  main  channel.  The  last  set 
of  experiments  were  conducted  in  a  meandering  channel,  in  a 
straight  valley. 


2  .  2  Exp  er imen t al  F 1 ume 

The  experiments  were  conducted  in  a  tilting  rectangular 
channel  4  feet  (1.22  m)  wide,  3  feet  (0.90  m)  deep,  and  60 
feet  (18.29  m)  long,  with  adjustable  slope.  The  side  walls 
of  the  flume  are  made  of  plexiglass.  The  plywood  bed  is 
covered  with  a  smooth  layer  of  fiberglass  resin.  This  flume 
was  used  as  a  base.  All  channels,  consructed  from  concrete, 
were  placed  within  this  rectangular  flume.  Considerable  care 
was  taken  to  ensure  that  the  concrete  blocks  were  uniformly 
laid  on  the  bed  of  the  flume  and  the  bed  was  level  prior  to 
adding  roughness  or  tilting  the  flume. 

The  flume  was  supported  on  three  sets  of  synchronized 
screws  which  could  be  driven  by  an  electric  motor  for  slope 
adjustment.  Water  was  pumped  from  a  sump  into  the  head-tank 
of  the  flume  through  a  12  inch  (30.46  cm)  overhead  pipe. 
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fitted  with  an  orifice  meter,  from  which  the  discharge  was 
obtained.  The  flume  was  fitted  with  a  tilting  tailgate  for 
control  of  the  depth  of  flow. 

2.2.1  Boundary  Roughnes 

For  most  of  the  testing,  the  finished  smooth  surface  of 
the  concrete  was  used.  In  some  tests,  half  inch  (12. 7  mm ) 
diameter  hemispheres,  closely  packed,  were  used  as  the 
roughness  elements.  This  roughness  was  fabricated  from  hard 
rubber  into  rectangular  sheets  using  a  steel  mold.  The 
roughness  height  is  0.0208  feet  (6.35  mm).  These  rubber  mats 
were  placed  in  the  flume  for  a  length  of  about  30  feet  (9.15 
m)  and  the  velocity  and  shear  measurements  were  carried  out 
at  a  section  located  at  least  20  feet  (6.1  m)  downstream 
from  the  start  of  the  roughness. 

2.2.2  Temp  e  r a  t  ur  e  Me  asurements 

Water  temperature  was  measured  in  order  to  determine 
the  viscosity.  The  temperature  was  taken  at  fixed  time 
intervals  during  tests  using  an  electronic  thermometer  which 
was  read  to  the  nearest  0.01  degree  Celsius.  Water 
temperature  ranged  from  17.0  degrees  to  19.0  degrees 
C  e  Is  iu  s  . 


2  .  3  Exp  er iment al  Channel 

All  compound  cross-sections  simulating  the  main  channel 
and  the  flood  plain  were  made  of  concrete  which  was  steel 
trowelled  to  ensure  smooth  surfaces. 
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2.3.1  Symme  t  r ical  Channel 

The  compound  cross-section  simulating  the  main  channel 
and  symmetrical  flood  plain  was  made  of  two  rows  of  concrete 
blocks.  For  the  first  seven  experiments,  the  depth  of  the 
main  channel  h  below  the  flood  plain  level  was  about  A. 3 
inches  (10.92  cm)  and  for  the  later  experiments  h  was 
reduced  to  3  inches  (7.62  cm)  by  placing  a  layer  of  concrete 
in  the  main  channel.  The  width  of  the  main  channel  B  was 
maintained  at  8  inches  (20.32  cm)  for  all  these  experiments. 
The  total  width  of  the  flood  plain  b  was  equal  to  40  inches 
(1.02  m)  in  all  these  experiments. 

The  measurements  were  made  in  the  central  part  of  the 
channel  at  a  distance  of  approximately  30  feet  (9.15  m. )  from 
the  entrance,  after  making  certain  that  the  flow  was  fully 
developed  and  that  the  depth  of  the  flow  was  approximately 
constant.  For  the  first  seven  experiments,  due  to  the 
presence  of  a  small  hump  in  the  bed  of  the  main  channel,  the 
effective  slope  for  the  main  channel  could  be  assessed  only 
approximately.  The  slope  given  in  Table  1  for  experiments 
1  -  7  is  to  be  taken  as  the  slope  of  the  flood  plain.  For 

all  further  experiments,  the  effect  of  this  hump  was  removed 
by  adding  a  layer  of  thick  cement  to  the  main  channel  bed 
and  by  making  it  parallel  to  the  flood  plain  bed. 

The  velocity  measurements  were  made  in  one  half  of  the 
compound  section.  Measurements  of  the  bed  shear  stress  in 
the  main  channel  and  flood  plain  were  made  by  using  the 
Prandtl  tube  as  a  Preston  tube  (Preston,  1954,  Rajaratnam, 
1965,  and  Rajaratnam  &  Muralidhar,  1968).  A  series  of  seven 
additional  experiments  were  done  mainly  to  obtain  more  data 


on  the  bed  shear  stress  distributions  in  the  main  channel 
and  flood  plain.  In  the  first  eight  experiments  (the 
significant  details  are  recorded  in  Table  1),  the  boundary 
surfaces  were  generally  smooth.  For  the  ninth  experiment, 
the  bed  of  the  flood  plain  was  made  rough  by  spreading  the 
rubber  mat  with  hemispherical  projections  of  0.5  inch  (12.7 
mm)  diameter. 

In  addition,  a  number  of  special  experiments  were 
performed  with  the  flow  confined  to  the  main  channel,  with 
smooth  and  rough  boundaries  to  obtain  some  supplementary 
information  on  the  distribution  of  the  bed  shear  stress  in 
rectangular  channels  of  small  and  medium  aspect  ratios. 

These  result  will  be  presented  in  Appendix  A. 

2.3.2  Asymmetrical  Channel 

In  order  to  asses  the  effect  of  aspect  ratio  in  the 
main  channel  and  flood  plain  channel  three  asymmetrical 
channels  were  tested.  In  addition,  a  multi-stage  flood  plain 
channel  was  also  constructed. 

2  .  3  .  2  .  1  As  ymme  t r ica 1  Channel  Type  A 

One  side  wall  of  the  main  channel  was  extended  by  a  row 
of  concrete  bricks.  The  depth  of  main  channel  h  below  the 
flood  plain  level  was  about  3  inches  (7.62  cm).  The  width  of 
the  main  channel  B  was  maintained  at  8  inches  (20.3  cm).  The 
flood  plain  width  b  was  equal  to  20  inches  (50.8  cm),  see 
Figure  (2.1A).  The  velocity  and  shear  measurements  were  made 
for  the  whole  cross-section  at  a  section,  half  way  down  the 
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2  »  3  .  2  .  2  As  ymme  t  r  ic  a  1  Ch  annel  Type  B_ 

Since,  in  the  previous  experiments,  the  width  of  the 
main  channel  was  not  large  enough  and  the  depth  of  the  main 
channel  from  the  flood  plain  level  h  was  high,  the  concrete 
blocks  were  removed  (See  Figure  2. IB).  The  depth  of  main 

channel  h,  was  1.25  inches  (3.2  cm), the  width  of  main 

channel  20  inches  (50.8  cm)  and  the  width  of  flood  plain  was 
8  inches  (20.3  cm).  Measurements  were  conducted  30  feet  from 
the  head-t  ank . 

2  .  3  .  2  .  3  Asymme  t  r  ic  al  Channel  Type  C_ 

In  order  to  simulate  reasonably  wide  main  and  flood 
plain  channels,  the  bottom  of  the  main  channel  was  filled  up 
with  mortar.  The  depth  of  main  channel  h  below  the  flood 
plain  level  was  about  3  inches  (7.6  cm).  The  width  of  the 
main  channel  B  was  maintained  at  28  inches  (71.1  cm).  The 

flood  plain  width  b  was  eaual  to  20  inches  (50.8  cm). 

Five  experiments  were  conducted.  The  middle  of  the 
flume  was  chosen  for  cross-sectional  measurements.  The 
velocity  and  shear  stress  measurements  were  made  over  the 
entire  compound  cross-section.  Figure  2. 1C  shows  the 
dimensions  of  the  channel  and  Table  2  shows  the  significant 
details  of  experiments  1  to  5. 

2  .  3  .  2 . 4  As  ymme  t  r  ica  1  Channel  Type  D_ 

In  addition,  one  test  with  multi-level  flood  plain  was 
conducted.  The  depth  of  main  channel  h,  was  1.25  inches  (3.2 
cm)  and  the  depth  of  first  level  of  flood  plain  h',  was  2.75 
inches  (7.0  cm).  The  main  channel  width  B  was  20  inches 


(50.8  cm).  The  first  level  flood  plain  width  was  b/2=8 
inches  (20.3  cm)  and  the  second  level  flood  plain  width  was 
b  '  /  2  =  2  0  inches  (50.8  cm).  Figure  2. ID  shows  the  detailed 
dimensions  of  this  channel. 


2.3.3  Meandering  Channel 

In  designing  a  meandering  channel  that  is  to  be 
representive  of  a  typical  meandering  river,  one  has  to 
consider  several  parameters.  The  parameters  involved  can  be 
classified  into  four  groups: 

1.  The  ch a nne 1 -g e ome t ry  characteristics, 

2.  the  flow  characteristics, 

3.  the  fluid  properties, 

4.  the  sediment  properties. 

Obviously,  these  are  too  many  variables  to  be  handled 
in  the  limited  size  of  laboratory  flume  described 
previously.  Therefore,  sediment  was  eliminated  from  the 
study;  clear  water  was  used  as  the  fluid;  the  movable  bed 
was  replaced  by  a  smooth  rigid  boundary;  the  channel  was 
built  with  a  uniform  rectangular  cross-section.  The 
meandering  channel  with  constant  curvature  was  constructed. 
Short  tangents  connected  the  reverse  curves  with  one 
another . 

With  this  simplified  model  (See  Figure  2.2),  the 
geometry  of  the  curved  channel  can  be  defined  by  the  central 
angle  of  the  bend  9C ,  centerline  radius  R,  the  tangent 
length  T,  the  width  of  the  channel  B,  and  the  depth  of  flow 

D  . 


B.C.  Yen  (1965)  conducted  a  survey  of  the  most 
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Figure  2.3  Location  of  Cross-Sections  in  Curved  Channel 
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important  experimental  work  on  flow  in  open  channel  bends. 
Based  on  this  survey,  he  maintained  that  those  channels  had 
a  few  sho r t -c omings .  These  s ho r t -c oming s  include  the  fact 
the  aspect  ratio  B/D  was  small,  or  if  the  channel  had  a 
large  width,  the  total  central  angle  ©c  is  too  large  or  too 
small.  However,  in  most  cases  there  was  only  one  bend. 
Leopold  &  Langbein  (1966)  used  a  statistical  analysis  of 
rivers  to  give  a  range  for  sinuosity,  or  tightness  of  bend, 
which  is  expressed  as  the  ratio  of  the  length  of  the  channel 
in  a  given  curve  lw,  to  the  wave  length  A,  between  1 .3:1 
and  4:1.  Another  property  is  the  ratio  of  wave  lengti  of  the 
curve  and  the  radius  of  curvature.  The  appearance  of 
regularity  depends  in  part  on  how  constant  this  ratio  is. 
This  value  may  change  from  the  water  surface  to  the  bottom 
of  the  channel.  The  average  value  of  1W/A  is  about  4.7:1. 
Also  Leopold,  Pagnold  &  Wolman  (1960)  gave  an  average  value 
of  R/B  of  2.3  for  many  rivers  and  channels,  and  an  average 
value  of  3.24  for  the  few  data  they  obtained  for  Mississippi 
River. 

Based  on  all  this  information,  and  the  limitation  of 
the  size  of  the  available  flume  in  the  laboratory,  the 
planimetric  geometry  of  the  model  was  chosen  with  9C=90 
degrees,  R/B=3.0,  A/R=3.02,  and  1W/A  =1.13.  The  model  was 
built  with  the  following  dimensions:  B=10  inches  (25.4  cm), 
R=30  inches  (76.2  cm),  tangential  length  T=4  inches  (10.2 
cm),  and  h,  the  depth  of  the  curved  channel  bed  below  the 
flood  plain  equal  to  1.5  inches  (3.R  cm).  Seven  identical  90 
degree  curves  of  reversed  direction  were  connected  by  the  4 
inches  (10.16  cm)  straight  reach,  where  it  formed  total  of 
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three  loops.  It  was  believed  that  developed  flow  would  be 
obtained  in  the  second  loop.  The  channel  was  60  feet  (18.3 
m)  long.  The  first  18  feet  (5.5  m)  was  a  straight 
symmetrical  main  channel  with  flood  plain  width  (len=18 
feet,  5.5  m)  of  19  inches  (48.3  cm)  on  either  side.  The 
middle  part  26.6  feet  (8.1  m)  of  the  flume  contained  three 
loops,  and  the  rest  was  straight  channel  with  a  tail  gate  at 
the  downsrteam  end  (lex=13.5  feet  or  4.1  m) .  The  general 
layout  of  the  model  is  shown  in  Figure  2.2.  As  shown  in  this 
schematic  plan  view,  the  amplitude  of  the  wave  was  20  inches 
(50.8  cm)  and  the  narrowest  width  of  the  flood  plain  in  the 
meandering  portion  was  9  i n ch e s ( c  = 2 2 . 9  cm). 

The  velocity  and  shear  measurements  were  made  in 
cross-sections  at  every  15  degree  increment  of  9  in  the 
second  loop  and  at  the  entrance  of  the  loop.  Figure  2.3 
shows  a  detailed  plan  view  of  these  sections.  The  vertical 
profiles  were  located  at  0.1  feet  (3.05  cm)  intervals  in  the 
flood  plain  and  at  1  inch  (2.54  cm)  intervals  in  the  main 
channel.  Velocities  were  measured  every  0.01  feet  (3  mm)  in 
vertical  profiles. 


2 . 4  Velocity  Me  a  sur  erne  n  t  s 

Prandtl  and  yaw  probes  were  used  for  measuring  the 
velocity  field.  In  straight  channels  with  very  weak 
secondary  currents,  the  velocity  measurements  were  made 
using  a  3  mm  or  1.5  mm  Prandtl  type  Pitot-static  tube.  The 
differential  pressure  was  measured  with  a  Pace  model  P1D 
variable  reluctance  differential  pressure  transducer  with 
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one  inch  of  water  for  full  scale  diaphragm  (zero  to  ten 
volts  on  the  demodulator  amplifier  respectfully).  Initially, 
readings  were  obtained  from  a  direct  reading  Pace  indicator 
(model  CD  25).  Later  on,  the  outputs  were  recorded  on  a 
strip  chart  recorder  and  averaging  was  done  by  hand. 

The  calibration  procedure  consisted  of  applying  static 
differential  heads  to  the  Prandtl  tube  and  adjusting  the 
demodulator  for  a  linear  output  of  zero  to  one  inch  of 
water.  According  to  the  manufacturer,  the  precision  of  these 
transducers  is  +0.1%. 

Each  modulator  was  provided  with  a  capacitor  which 
damped  the  circuits  to  reduce  the  frequency  responses  of  the 
output  to  a  level  suitable  for  input  to  the  strip  chart 
recorder.  Instantaneous  differential  pressure  between 
dynamic  and  static  heads  was  recorded  for  one  minute. 
Whenever  the  depth  of  flow  on  the  flood  plain  was  small,  the 
requirement  of  depth  of  flow  being  at  least  five  times  the 
outside  diameter  of  the  probe  could  not  be  met  with  a  3mm  CD 
probe.  Therefore,  for  small  depths  of  flow,  the  1.5  mm  OD 
Prandtl  Pitot-static  tube  was  used. 


2 . 5  Veloc i ty  Vector  Measurement 

For  velocity  measurements  in  the  curved  channel  with  a 
flood  plain,  the  conventional  Pitot  tube  would  not  suffice 
as  the  angle  of  attack  would  be  very  large.  In  addition,  the 
present  investigation  was  aimed  at  measuring  all  the  three 
components  of  a  velocity  vector  which  would  require  the 
measurement  of  the  yaw  and  pitch  angles.  The  obvious  choice 
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was  a  five  hole  probe.  After  a  preliminary  study  in  the 
curved  channel  it  was  found  that  the  vertical  component  of 
the  velocity  vector  was  rather  small  compared  to  the 
horizontal  component  of  the  velocity  vector  and  therefore 
could  be  neglected.  Besides,  the  outside  diameter  of  five 
hole  probe  was  large  and  could  not  meet  the  requirement  of 
the  depth  of  flow  being  at  least  five  times  the  probe 
diameter.  Further,  the  complexity  of  the  calibration  curves 
ruled  out  the  possibility  of  using  the  five  hole  probe.  In 
two  dimensional  flows,  when  the  velocity  vector  is  in  the 
horizontal  plane  of  the  probe,  its  magnitude  and  direction 
can  be  conveniently  determined  by  means  of  the  yaw  probe. 

The  yaw  probe  was  prepared  from  stainless  steel  tubes 
of  1.5  mm  outside  diameter  and  had  a  45  degree  chamfer.  At 
the  beginning  of  the  experiment  a  rotating  mechanism  was 
attached  to  the  yaw  probe.  By  rotating  the  yaw  probe  until 
the  differential  pressures  between  center  tube  and  the  side 
tube  became  the  same,  it  was  possible  to  obtain  the 
direction  and  magnitude  of  the  velocity  vector.  The  lag  time 
between  rotating  the  yaw  probe  and  the  response  of  the 
differential  pressures  kept  the  yaw  probe  in  constant 
motion.  The  fluctuation  of  direction  of  the  velocity  vector 
added  a  greater  problem.  In  order  to  overcome  these 
difficulties,  an  alternative  method  was  used.  Ra jaratman  & 
Muralidhar  (1968)  gave  the  necessary  calibration  curves  for 
computation  of  magnitude  and  direction  of  the  velocity 
vector.  The  principle  involved  in  this  method  of  velocity 
measurement  is  described  briefly  in  the  following 
p  aragr aphs . 
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For  two-dimensional  flows,  if  at  any  point  the  velocity 
is  V  at  any  angle  of  attack  of  0  and  the  static  piezometric 
head  is  hQ,  the  piezometric  heads  indicated  by  the  three 
tubes  of  the  probe  set  at  that  point  are  written  as: 
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+  k 


(2.1) 
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and 


+  k  vf 

3  2g 


(2.3) 


where  k^,  k^,  and  k^  are  the  calibration  coefficients  and 
are  functions  of  0  mainly,  for  a  given  probe  neglecting 
viscous  effects  and  other  minor  correction  factors  (see 
Figure  2.4).  The  mean  calibration  curves  for  k^,  k^  and  k^ 
were  obtained  through  a  series  of  experiments,  with  0 
ranging  from  -60  to  60  degrees. 

Using  the  three  calibration  factors  k  ^  ,  k^  an  k  ^  a 
fourth  factor  k  is  defined  as: 


k  = 


(2.4) 


and  it  could  be  seen  that  k  is  mainly  a  function  of  0.  Using 
the  mean  curves  presented  in  Figure  2.4,  the  variation  of  k 
with  0  is  shown  in  Figure  2.5. 

Equation  2.4  in  combination  with  equations  2.1  to  2.3 
could  be  rewritten  as: 
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(2.5) 


From  equation  2.5,  k  is  found  and  hence  the  angle  of  the  yaw 
is  found  by  referring  to  Figure  2.5.  It  could  be  shown  that 
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(2.6) 


wh  ere 


(2.7) 


and 

Ah  =  h  -  h2  .  (2.8) 

It  could  be  shown  that  k.  is  only  a  function  of  9,  and  its 

variation  with  9  is  shown  in  Figure  2.6.  Using  this  figure, 

k  could  be  found  for  the  particular  angle  of  attack  9  and 
o 

the  magnitude  of  the  velocity  vector  is  determined. 


2 . 6  Boundary  Shear  Stress  Measurements 

The  Preston  tube  is  a  total  head  tube  resting  on  a 
boundary  and  facing  the  flow,  with  the  tube  in  a  zone  where 
the  law  of  the  wall  is  valid.  This  requirement  was  met  at 
all  times  and  whenever  the  depth  of  flow  was  shallow,  the 
1.5  mm  Pitot-static  probe  used.  The  Preston  technique  is  a 
convenient  method  of  measuring  the  boundary  shear  stress. 
The  differential  ^p  pressure  between  the  dynamic  and  the 
static  pressures  is  correlated  with  the  boundary  shear 
stress  T  .  These  correlation  was  developed  by  Preston 
(1954),  mainly  from  similarity  considerations  and  the 
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existence  of  the  law  of  the  wall.  This  correlation  is  given 
in  the  form  of 


AP  d‘ 

4  p  V 
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4  p  v'1 
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(2.9) 


After  Preston  (1954)  a  number  of  investigators  conducted 
several  experiments  to  check  Preston's  original  calibration, 
and  to  extend  the  use  of  the  technique  to  developing 
boundary  layers  and  to  study  the  effect  of  pressure 
gradients  (Hsu,  1955,  and  Bradshaw  &  Gregory,  1958). 
Pajaratnam  (1965)  defined  calibration  curves  theoretically 
by  assuming  a  suitable  velocity  distribution  in  the  wall 
region.  The  numerous  smooth  boundary  calibrations  showed 
only  a  few  percent  difference.  Rajaratnam  &  Furalidhar 
(1968)  showed  that  the  Prandtl  tube  could  be  used  as  a 
Preston  tube.  In  the  present  investigation,  Patel's  curves 
(Patel,  1965)  were  used  to  convert  ^p  to  TQ  and  the  required 
calibration  curves  for  the  3  mm  and  1.5  mm  Pitot-static 


probe  were  prepared 


2.6.1  Y  aw  Probe  a  s  a_  Preston  Tube 

In  the  meandering  channel,  the  direction  of  the  bed 
shear  stress  is  not  known  beforehand.  In  principle,  its 
magnitude  and  direction  can  still  be  found  by  rotating  the 
Preston  tube  till  the  maximum  dynamic  pressure  is  indicated. 
This  method  of  rotating  the  probe  is  time  consuming.  In 
order  to  overcome  these  difficulties,  Rajaratnam  and 
Muralidhar  (1968)  used  the  yaw  probe  on  the  boundary  and 
gave  the  necessary  calibration  curves  for  obtaining  the 


■ 


. 


direction  and  magnitude  of  the  bed  shear  stress.  The  same 
yaw  probe,  described  previously  in  connection  with  velocity 
measurement,  was  used  as  the  Preston  tube. 

Combining  various  expressions  in  equation  2.4  and  2.5 
the  following  can  be  written: 


P1  -  Po  +  klf 

p2  =  p0  +  k2( 

p3  =  p0  +  k3( 
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and  p^,  p^  and  p^  are  respectively  the  pressure  at  hole  1,  2 

and  3  . 

By  a  suitable  combination,  the  expressions  in  2.5  can 
be  written  as  the  following: 
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(2.14) 


In  measuring  shear  stress  at  any  point,  p^,  p^  and  p^  are 
known  quantities  so  that  k  could  be  computed  and  then  9, 
could  be  found.  Also  k^  and  k2  could  be  found  from  Figure 
2.4.  Since  all  the  quantities  on  the  left  hand  side  of 
equation  2.14  are  known,  f  (  T  +)  can  be  determined.  Using 
Patel's  calibration  curves  Tq  +  and  hence  tq  were  computed. 
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2 . 7  Depth  and  Slope  o  f  Water  Surface 

Depth  was  determined  from  manual  gauge  readings  of  the 
water  surface  along  the  flume  centerline  at  every 
cross-section.  These  readings  were  referenced  to  their  bed 
level  datum.  For  subcritical  flow  conditions,  the  tailgate 
was  adjusted  to  establish  uniform  flow  by  a  trial  and  error 
process.  The  bed  slope  was  set  and  uniform  flow  conditions 
were  approached  by  manipulating  the  tailgate  while  observing 
the  depth  of  flow.  They  were  measured  from  5  feet  to  50  feet 
Cl. 52  to  15.25  m)  from  head-tank  at  5  feet  intervals. 

After  a  few  runs  were  conducted,  it  was  found  that  the 
water  surface  measurement  with  a  point  gauge  was  not 
accurate  enough  for  this  study.  Seven  static  probes  were 
mounted  over  a  reach  betwen  15  to  45  feet  from  the  entrance 
at  intervals  of  5  feet.  The  differential  pressures  between 
the  static  probes  were  measured  with  a  pressure  transducer 
Pace  model  P1D.  Therefore,  the  uniform  flow  was  established 
by  manipulating  the  tailgate  until  the  6  differential 
pressures  between  the  static  probes  became  equal  and  the 
slope  could  be  calculated  by  dividing  the  differential 
pressure  reading  over  the  distance  between  the  two  probes. 
Since  the  drop  of  the  water  surface  could  be  read  within  a 
precision  of  1/10, 000th  of  a  foot,  the  slope  of  the  water 
surface  could  be  measured  accurately.  In  all  cases,  care  was 
taken  to  see  that  the  bed  and  the  water  surface  slopes  were 
essentially  the  same. 
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2  •  8  Special  Instrumentation 


In  the  straight  channel  a  y-z  electric  coordinator  and 
a  one  channel  strip  chart  recorder  were  quite  sufficient  to 
take  data  manually.  But,  in  the  case  of  the  meandering 
channel,  obtaining  data  manually  was  quite  time  consuming. 

In  addition  to  that,  processing  the  data  without  computer 
help  and  plotting  or  mapping  isovels  was  a  great  task  by 
hand. 

2.8.1  Coordinate  Pos  it  ions  and  Input  Panel 

The  yaw  probe  was  positioned  within  the  flume  by  a 
x-y-z-9  electric  coordinate  system,  where  9  is  the  degree  of 
rotation  of  the  yaw  probe  with  respect  to  the  x  direction. 

The  origin  of  the  x-y-z  coordinate  was  selected  in  such  a 
way  that  x=0.0  at  cross-section  NO.  1,  z=C.O  at  the  right 
hand  side  corner  of  the  flood  plain  and  y  or  y'=0.0  upward 
from  the  bed.  The  location  and  rotation  of  the  yaw  probe  (x-y 
-z-9)  could  be  controlled  from  a  remote  control  switch 
box,  which  in  turn  drove  a  small  DC  electric  motor  to  a 
desired  position.  The  coordinate  voltages  were  zeroed  with 
the  reference  to  the  flume  origin.  The  DC  voltages 
representing  the  coordinates  were  then  fed  into  the  data 
acquisition  system. 

This  system  also  was  controlled  from  a  small  panel 
where  coordinate  datum  planes  could  be  set  and  calibrations 
performed.  This  equipment  is  indicated  in  the  photograph  in 


Figure  2.7. 
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Data  Aquisition  System 


Meandering  Channel 


Yaw  Probe 


Figure  2.7  Photograph  of  Data  Aquisition  System 
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2.8.2  Data  Acquisition  System 

The  acquisition  system  was  made  up  of  the  following 
elements,  which  are  represented  in  Figure  2.8. 

1.  20  Channel  Scanner  (Vidar  604) 

2.  Integrating  Digital  Voltm  (Vidar  500) 

3.  Data  Coupler  (Control  Equipment  Corporation  model  310) 

4.  Scan  Counter 

5.  Magnetic  Tape  Recorder  (Techtran) 

All  inputs  for  the  data  acquisition  system  were  in  the  range 
of  -20.0  to  20.0  volts  and  each  was  assigned  to  a  particular 
channel  of  the  Scanner.  The  Scanner  sequentially  sampled 
selected  channels,  routing  the  analogue  signals  to  the 
Digital  Voltm  where  they  were  digitized  and  converted  to  a 
Binary  Coded  Decimal  (BCD)  Code.  The  Data  Coupler  formulated 
the  BCD  signal  and  routed  it  to  the  magnetic  tape  recorder 
(Techtran).  The  Techtran,  on  write  mode,  stored  the 
information  on  a  casset  tape  and,  at  a  later  period,  the 
tape  could  be  read  to  PDP8  computer. 

The  desired  number  of  scan  cycles  were  selected  on  the 
Scan  Counter  and  the  scan  rate  was  selected  on  the  Scanner. 
After  several  trials,  2400  characters  per  minute  gave  the 
best  result,  with  approximately  80  samplings  per  minute  in 
each  channel . 

Through  a  Basic  language  program,  the  information  on 
the  magnetic  tape  was  read  to  the  Digital  Computer  (PDP8) 
and  the  data  was  stored  on  the  Flupy.  Furthermore,  PDP8 
through  several  programs,  was  able  to  convert  each  set  of 
readings  to  the  desired  average  of  that  reading  and,  if 
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Figure  2.8  Flow  Chart  of  Equipments 
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necessary,  by  going  through  a  subroutine  and  applying 
calibration  curves  to  obtain  the  value  of  the  variables. 
Then,  this  information  from  the  Hydraulic  Laboratory  was 
transferred  to  the  University  Computer  Center.  From  there 
the  data  were  stored  in  suitable  matrices.  It  could  be 
called  for  under  the  Michigan  Time  Sharing  System  (MTS) 
using  an  Amdahl  4  7  0  V  /  6  .  All  data  processing  and  analyses 
were  carried  out  using  Amdahl  4  7  0  V  /  6  and  the  results  were 
plotted  using  a  new  version  of  Calcomp  plotter.  Surface  II 
Program  (Sampson,  1978)  had  more  capability  for  mapping, 
therefore  it  was  decided  to  use  this  program  for  i  s  o  v  1 
ma  p  s  . 

2.8.3  Test  Method 

The  following  procedure  was  normally  carried  out  to 
complete  a  single  test: 

1.  The  channel  slope  was  selected, 

2.  the  flume  slope  was  set, 

3.  a  discharge  was  selected, 

4.  uniform  flow  was  established  by  measuring  water  depth, 
and  adjusting  the  tailgate  elevation  until  the 
differential  pressure  between  each  pair  of  the  static 
probes  became  equal, 

5.  final  depths  were  measured,  the  water  surface  was 
plotted,  and  slope  was  determined, 

6.  the  velocity  distribution  was  measured  by: 

a.  selecting  the  cross-section, 

b.  zeroing  and  calibrating  the  coordinate  positioner, 
and  the  pressure  transducers,  and  changing  the  angle 
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of  the  yaw  probe  if  necessary, 

7.  making  sure  scan  rate  is  set  at  80  cycles, 

8.  positioning  the  velocity  probe,  selecting  the  channels 
to  be  scanned,  and  starting  auto  counter, 

9.  recording  the  data  on  magnetic  tape, 

10.  shear  measurements  were  carried  out  by  following  the 
same  procedure  as  above  with  the  yaw  probe  positioned  on 
the  bounda  ry  . 

2  .  9  F 1 ow  Visualisation 

Dye  was  injected  parallel  to  the  flow,  by  means  of  a  1 
mm  OD  tube  connected  to  a  small  constant  pressure  head-tank. 
Dispersion  of  the  dye  was  photographed  from  a  fixed  position 
approximately  4  feet  (1.2  m)  above  the  flow.  These  photos 
served  many  purposes  such  as  measuring  the  rate  of  expansion 
of  the  width  of  the  dye  plume,  or  the  extent  of  the 
intrusion  of  the  main  channel  flow  into  the  flood  plain 
region. 


. 


3.  Interaction  Between  Straight  Main  Channel  And  Flood  Plain 

Flows  in  a  Symmetrical  Channel 
A  preliminary  study  was  done  in  a  narrow  main  channel  with 
wide  symmetrical  flood  plain. 

The  construction  and  dimensions  of  this  symmetrical 
flume  were  described  in  detail  in  Chapter  II,  Section  2.3.1. 
Originally,  it  was  used  to  check  the  basic  idea  of  exchange 
of  momentum  between  the  main  channel  and  the  flood  plain. 


3  .  1  Concep  t  o  f  Momentum  Exchange  and  Interaction  P  e  g i on 

In  a  compound  cross-section  channel,  composed  of  two 
similar  channels  with  the  same  longitudinal  slope,  the  one 
with  a  larger  depth  of  flow  at  any  given  level  will  have  a 
higher  velocity  than  the  one  with  a  lower  depth  of  flow  at 
the  same  level.  Therefore,  it  is  obvious  in  a  compound 
cross-section  like  a  river  with  a  flood  plain,  at  any  level 
above  the  flood  plain,  the  main  channel  will  have  higher 
velocity  because  of  greater  depth  of  flow.  It  can  be  seen  at 
any  of  these  levels  that  there  exists  a  lateral  velocity 
gradient,  higher  velocity  in  the  middle  of  the  main  channel 
and  lower  in  the  flood  plain.  The  result  of  this  difference 
will  be  that  a  part  of  the  main  flow  near  the  flood  plain 
will  slow  down  and  a  part  of  the  flow  in  the  flood  plain 
near  the  main  channel  will  speed  up.  Perhaps  this  action 
will  occur  in  such  a  way  that  the  velocity  of  the  main  flow 
decreases  smoothly  toward  the  flood  plain  until  it  reaches 
the  flood  flow  velocity.  This  is  the  result  of  transport  of 
the  longitudinal  momentum  to  the  flood  plain  and  the 
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vertical  plane  between  the  main  channel  and  flood  plain  will 
experience  appreciable  turbulent  shear  stress  in  the 
longitudinal  direction. 

In  order  to  test  the  above  ideas,  data  were  taken  from 
the  isovel  contours  given  by  Sellin  (1964)  and  plotted.  A 
very  simple  calculation  resulted  in  undisturbed  velocity  for 
the  main  channel  and  the  flood  plain  channel.  Comparison  of 
these  velocities  with  those  obtained  from  the  isovel  showed 
main  channel  velocities  were  reduced  and  flood  plain 
velocities  increased. 

In  order  to  be  able  to  pursue  the  forgoing  idea,  the 
symmetrical  channel  was  tested  for  different  depths  of  main 
channel  and  flood  plain  flows. 


3  .  2  Experiments  An  d  Experimental  P  e  s  u 1 1  s 

Measurements  were  made  in  the  central  part  of  the 
channel  (at  a  distance  of  approximately  30  feet  (9.15  m) 
from  the  entrance)  after  making  certain  that  the  flow  was 
fully  developed  and  that  the  depth  of  flow  was  approximately 
constant.  As  previously  explained,  for  the  first  seven 
experiments,  the  effective  slope  for  the  main  channel,  could 
only  be  approximated.  The  slope  given  in  Table  1  for 
experiments  1  -  7  is  to  be  taken  as  the  slope  of  the  flood 
plain. 

The  velocity  and  shear  stress  measurements  were  made  in 
one  half  of  the  compound  section.  A  series  of  seven 
additional  experiments  were  made,  mainly  to  obtain  more  data 
on  the  bed  shear  in  the  main  channel  and  the  flood  plain.  In 
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TABLE  1 

Significant  Details  of  the  Experiments 


B  =  0.67  feet,  b  =  3.33  feet 


Exp  t . 

h  (f  t ) 

D(f  t) 

d  (f  t) 

0 

S 

D/d 

F 

Shear 

No  . 

Ccf  s  ) 

x  1  0  3 

lb  s  x 1 0  3 

sq.ft 

1 

0.346 

0.483 

0.137 

0.98 

0.70 

3.53 

0.52 

18.37 

2 

0.358 

0.4  10 

0.052 

0.50 

0.40 

7.88 

0.5  9 

90.96 

3 

0.358 

0.584 

0.226 

2.00 

1.14 

2.58 

0.59 

16.56 

4 

0.358 

0.538 

0.180 

1.50 

1.07 

2.99 

0.56 

17.59 

5 

0.358 

0.407 

0.12° 

0.58 

0.29 

3.16 

0.3  1 

23.32 

6 

0 .358 

0.518 

0.160 

0.18 

0.03 

3.24 

0.08 

1.50 

7 

0.358 

0.450 

0.092 

0.18 

0.04 

4.89 

0.13 

8.50 

8 

0.245 

0.360 

0.115 

0.54 

1.18 

3.13 

0.39 

2.79 

9  * **  * 

0.271 

0.338 

0.067 

0.32 

0.94 

5.04 

0.22 

- 

10  * 

0.245 

0.550 

0.305 

— 

0.56 

1.80 

— 

2.07 

1  1  * 

0.245 

0.498 

0.253 

— 

0.53 

1  .  9  7 

— 

2.08 

1  2  * 

0.245 

0.422 

0.177 

— 

0.62 

2.38 

— 

1.32 

1  3  * 

0.245 

0.375 

0.130 

— 

0.70 

2.88 

— 

3.2  7 

14  * 

0.245 

0.345 

0.100 

— 

0.82 

3.45 

— 

3.67 

15  * 

0.245 

0.328 

0.083 

— 

0.83 

3.95 

— 

5.02 

16  * 

0.245 

0.264 

0.019 

1.27 

13.89 

—  — 

38.42 

*  Additional  experiments. 

**  Rough  flood  plain  0.5  inch  diameter  h emi -s ph er i ca 1 


roughness 


the  first  eight  experiments  (see  Table  1  for  the  significant 
details  of  the  experiments),  the  boundary  surfaces  were 
generally  smooth.  For  the  ninth  experiment  only,  the  bed  of 
the  flood  plain  was  made  rough.  In  addition,  a  number  of 
special  experiments  were  made  with  the  flow  confined  to  the 
main  channel,  to  obtain  supplementary  information  on  the 
distribution  of  the  bed  shear  in  rectangular  channels  of 
small  aspect  ratio.  These  results  are  presented  in  Appendix 
A. 

The  velocity  profiles  for  three  typical  experiments  are 
shown  in  Figure  3.2.  Figure  3.2(a)  shows  the  profiles  for  a 
large  Froude  number  F  of  0.52,  wherein  F  is  taken  as  the 
ratio  of  the  mean  velocity  to  the  square  root  of  the  product 
of  acceleration  due  to  gravity  and  the  depth  of  flow  D. 
Figure  3.2(b)  shows  the  profiles  for  F=0.08.  Figure  3.2(c) 
shows  the  results  for  the  ninth  experiment  which  had  a  rough 
bed  for  the  flood  plain.  Considering  the  flood  plain,  in 
Figure  3.2(a),  the  velocities  increase  continuously  as  one 
moves  towards  the  main  channel.  This  increase  can  be  noticed 
in  Figure  3.2(c)  and  to  a  lesser  extent  in  Figure  3.2(b).  In 
the  main  channel.  Figure  3.2(a),  the  velocity  increases  with 
y,  the  distance  from  the  bed  up  to  some  level,  and  decreases 
continuously  to  eventually  become  approximately  constant 
near  the  water  surface.  This  behaviour  is  noticeable  in 
Figure  3.2(b)  as  well  as  Figure  3.2(c). 

Figure  3.3(a)  shows  the  variation  of  u  with  z  from  the 
centerplane  for  different  horizontal  planes  (i.e.  different 
values  of  y).  From  Figure  3.3(a)  it  can  be  seen  that  u 
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(a)  (half)  Sectional  View 


(d)  Concept  of 

Flow  Reduction 


Figure  3.1  Definition  Sketches 


Expt  9  (d  =  0.067  ft) 
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Figure  3.2a,  b  &  c  Typical  Velocity  Profiles  (in  y  direction) 
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Figure  3.3a,  b  &  c  Typical  Velocity  Profiles  (in  z  direction) 
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decreases  from  a  maximum  value  of  u  on  the  centerline 

m 

rather  rapidly  to  a  value  of  u  ^  at  the  junction  plane 

between  the  main  channel  and  flood  plain  and  then  decreases 

continuously  from  u^  to  eventually  approach  the  undisturbed 

flood  plain  velocity  of  u.  Figures  3.3(b)  and  3.3(c)  show 

two  other  typical  profiles,  one  for  low  F  and  the  other  for 

a  rough  flood  plain  and  similar  variation  can  be  noticed.  In 

the  immediate  vicinity  of  the  junction  plane,  in  some 

experiments,  a  sharp  increase  in  velocity  was  noticed  (see 

for  example  Figure  3.3(b)).  This  sudden  increase  has  also 

been  reported  by  other  observers  (Goncharov,  1964). 

The  variation  of  the  bed  shear  stress  in  the  main 

channel  and  flood  plain  is  shown  in  Figure  3.4  for  the  eight 

experiments  as  well  as  the  additional  series.  In  the  main 

channel  the  bed  shear  decreases  from  a  maximum  value  of  t  ™ 

o  m 

to  zero  at  the  corner  of  the  banks  .  In  the  flood  plain, 
decreases  continuously  from  a  maximum  at  the  junction  plane 
to  approach  an  undisturbed  value  of  .  T0oo  agaan 

decrease  as  the  bank  of  the  flood  plain  is  approached. 


3 . 3  Analysis  o  f  Fxperi mental  Results 

Some  preliminary  calculations  indicated  that  in  the 

part  of  the  flow  above  the  bed  of  the  flood  plain  (i.e. 

y>h),  the  velocity  profiles  in  the  transverse  direction 

(i.e.  u(z))  if  viewed  with  respect  to  u  ,  the  undisturbed 

oo 

velocity  in  the  flood  plain,  appeared  to  possess  the 
property  of  similarity.  After  further  analysis,  it  was  found 
that  the  u(z)  profiles  could  be  split  into  two  parts,  one  in 
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Figure  3.4a  &  b  Bed  Shear  Stress  Profiles 


48 


the  main  channel  and  the  other  in  the  flood  plain.  If  u  and 

m 

are  the  values  of  the  longitudinal  velocity  in  the 
centerline  of  the  main  channel  and  at  the  junction  of  the 


main  channel  and  flood  plain  for  any  given  horizontal  plane 


for  which  u^  is  the  corresponding  undisturbed  flood  plain 


velocity,  it  was  found  that  if  (u  -  u  '  )  /  (  u  -  u')  is 

m  m  m 

plotted  against  n=z/bm  where  b^  is  the  length  scale,  taken 

as  equal  to  z  where  (u  -  u')/(u  -  u')=0.25,  the  profiles 

m  m  m 

are  approximately  similar  as  shown  for  a  few  typical 


experiments  in  Figure  3.5.  The  data  are  described  by  a  curve 


of  the  eq  ua  t io  n : 


f(n)  =  — - -2-,  =  1  -  o.75n2  (3.1) 

u  -  u 
m  m 


In  the  above  description,  u  and  u  '  are  the  reference 

m  m 

velocities  and  b  is  the  length  scale.  If  the  main  channel 

m 

was  sufficiently  wide  (Pajaratnam  &  Muralidhar  1969),  u^ 
could  be  represented  by  the  equation: 


u 

m 


y  u  . 

5.75  1  o  g  (  -  )  +  5.5 

V 


(3.2) 


where  u  =/~~g  D  S~ ,  S  is  the  slope  of  the  channel  and  V  is 
*  °  o  o 

the  kinematic  viscosity  of  water.  The  undisturbed  flood 
plain  velocity  u^  would  be  given  by  the  equation: 


u 


oo 


u 


=  5.75  log ( 


y  u* 


oo 


*00 


V 


)  +  5.5 


(3.3) 


where  u.  =  /g d  S  ,  and  y'  is  flood  plain  (see  Figure 


*oo  ~  o 


3.1) 
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y'  = 


*  0.063  ft  (Expt  9) 

*  0.070  (Expt  7) 

*  0.13  (Expt  6) 

*  0.12  (Expt  5) 

*  0.17  (Expt  4) 


0.4  0.6  0.8 

n  ~  z/bm 

Figure  3.5  Velocity  Profile  Similarity:  Upper  Part  of  Main  Channel 
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Then  it  could  be  shown  that: 


u  -  u 
m_ go 


u 
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(3.4b) 


wherein  C  is  a  constant  equal  to  9.06. 

In  this  study,  the  width  of  the  main  channel  is  not 
large  and,  as  a  result,  the  effects  of  mixing  at  the 
interface  of  the  main  channel  and  flood  plain  have  permeated 
to  the  centerline  of  the  main  channel.  Using  equation  3.4  as 
a  guide,  a  scheme  was  found  for  correlating  u^  by  plotting 

^um  “  Uoo)/(u00)  agai-nst  y ' /d  (see  Figure  3.6(a)).  For  any 

experiment,  (u^  -  Uoo^  ^Uoo^  decreases  with  y'/d  to 

eventually  approach  an  asymptotic  value  as  y'/d  approaches 

unity.  Using  the  results  in  Figure  3.6(a),  Figure  3.6(b)  was 

prepared  which  shows  that  (um  -  u00)/(uco)  for  any  given 

value  of  y'/d  increases  linearly  with  D/d.  In  any  flood 

plain  problem,  if  the  flood  plain  is  sufficiently  wide,  as 

was  (approximately)  the  case  in  the  present  experiments, 

then  u^  could  be  found.  The  behaviour  of  the  length  scale 

b^  is  studied  in  Figure  3.7(a)  wherein  b> ^ / d  is  plotted 

against  y'/d  for  different  values  of  D/d.  For  any  given 

value  of  D/d,  variation  of  b  /d  with  y'/d  is  not  large  and 

m 

if  an  average  value  is  adopted,  this  average  value  of 
b  /D  is  seen  to  increase  linearly  with  D/d  as  shown  in 

ID 

Figure  3.7(b)  and  is  described  by  the  equation: 
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Figure  3.7  Length  Scale:  Upper  Part  of  Main  Channel 
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The  velocity  profiles  in  the  flood  plain,  when  plotted 
with  (uffl  -  uQ0)/(u^  -  u^)  against  r,'  =  (z'/bf)  where 
z'  =  (z  -  B/2)  and  bf=z'  where  (um  -  u00)/(u^  -  u^^O.5,  were 
found  to  be  (approximately)  similar  and  a  few  typical  plots 
are  shown  in  Figure  3.8.  These  curves  are  satisfactorily 
described  by  the  error  function.  In  some  plots  of  Figure 
3.8,  there  is  appreciable  scatter,  particularly  for  the 
lower  most  level,  i.e.  for  the  Pitot  tube  sitting  on  the  bed 
of  the  flood  plain.  This  observation  will  predict  scatter  in 
bed  shear  profiles.  The  behaviour  of  the  velocity  scale 

^  u  m  “  uoo)  /  (u  )  is  studied  in  Figure  3.9(a)  wherein  it  is 
found  that  for  y'/d=0.3,  it  varies  mainly  with  D/d  and  this 
variation  is  shown  separately  in  Figure  3.9(b).  In  Figure 

3.9(b),  (um  ~  Uoo^^Uoo^  s  seen  to  increase  linearly  with 
D/d.  Figure  3.9(c)  shows  that  the  length  scale  b^/d  is 
(almost)  independent  of  y'/d  and  varies  mainly  with  D/d. 

This  variation  is  seen  to  be  linear  (see  Figure  3.7(b))  and 
could  be  described  by  the  equation: 


2.02  — 

d 


-  2.71 


(3.6) 


The  flow  in  the  main  channel  below  the  bed  level  of  the 
flood  plain  is  analysed  in  an  empirical  manner.  When  the 
velocity  profiles  in  this  region  were  plotted  with  u/uq 

against  y/h  where  uq  is  the  value  of  u  for  y=h,  as  shown  in 
Figure  3.10,  it  could  be  seen  that  for  approximate 
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Figure  3.9  Velocity  and  Length  Scales:  Flood  Plain 


Figure  3.10  Velocity  Profile  Similarity:  Lower  Part  of  Main  Channel 
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calculation  purposes,  an  average  curve  (shown  by  a 
continuous  line  in  Figure  3.10)  could  be  used.  For  this  mean 
curve,  the  maximum  value  of  u/uo=l.06  and  occurs  at 
y/h=0.5.  To  use  this  curve,  u^  could  be  obtained  from  Figure 
3.6  for  the  lowest  (experimental)  value  of  y'/d.  An 
alternate  method  of  analysing  the  flow  in  this  region  is 
indicated  later. 

From  the  plot  of  the  shear  stress  measurements,  in  the 
flood  plain,  (Figure  3.M  it  is  easy  to  see  that  the  bed 
shear  stress  in  the  flood  plain  increases  as  the  junction 
plane  with  the  main  channel  is  approached.  If  is  the 

undisturbed  shear  stress  in  the  flood  plain  (occurring  at 
points  sufficiently  away  from  the  banks  of  the  flood  plain) 
and  if  Tq  is  the  shear  stress  at  any  distance  z'  the  maximum 
value  of  At  =  (t«  -  T~  ),  written  as  At  occurs  at  the 

edge  of  the  flood  plain  (i.e.  for  z'=0.0).  A  plot  of  the 
variation  of  At  /  At  „  with  z'  in  terms  of  a  suitable  length 
scale  showed  that  for  very  approximate  purposes,  the 
distribution  of  this  dimensionless  excess  bed  shear  stress 
could  be  considered  to  be  similar  but  there  was  considerable 
scatter.  If  1.  is  the  value  of  z'  where  At  =0.0,  the 
variation  of  1^/d  with  D/d  is  shown  in  Figure  3.11(a).  From 
Figure  3.11(a),  it  can  be  said  that  l*/d  increases  almost 
linearly  with  D/d.  Figure  3.11(b)  shows  a  plot  of  the 
variation  of  T0m/T0  with  D/d.  Even  though  it  is  difficult  to 
draw  a  trend-line  in  Figure  3.11(b),  it  can  be  seen  that  for 
D/d  =  10.0,  Tom=2  T  .  Similar  large  increases  have  been 
reported  by  Myers  and  Elsawy  (1975). 

Considering  the  bed  shear  stress  in  the  main  channel. 
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Figure  3.11  Bed  Shear  for  Flood  Plain  and  Defect  in  Main  Channel 
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the  results  of  the  additional  experiments  as  well  as  those 
of  experiment  8  are  shown  in  Figure  3.11(a)  wherein  Tom  is 
the  corresponding  TQ  which  will  exist  if  the  side  walls  of 
the  main  channel  were  extended  to  a  total  depth  of  D.  A 
curve  for  the  variation  of  Tom/ p  g  D  Sq  with  B/D  based  on 
the  available  results,  with  additional  experiments  performed 
by  the  author,  is  given  in  Figure  3.12.  From  Figure  3.11(a) 
it  appears  that  the  maximum  bed  shear  stress  in  the  main 
channel  is  reduced  by  about  40%  because  of  the  interaction 
with  the  flood  plain  flow.  Similar  decreases  have  been 
reported  by  Myers  and  Elsawy  (1975).  It  was  also  found  (see 
Figure  3.11(c))  that  the  variation  of  T  (z)/t  m  in  the  main 
channel  with  z/B/2  was  approximately  the  same  as  that  of 
rectangular  channels  without  flood  plains.  Fence,  the  main 
effect  appears  to  be  a  reduction  in  the  magnitude  of  the  bed 
shear  in  the  main  channel. 

If  is  the  average  shear  stress  exerted  by  the  main 

channel  flow  on  the  flood  plain  flow  at  the  junction  plane, 

(T*  d)  could  be  recognized  as  equal  to  (Integral  of  At  q  d  z  ' 
from  0  to  z' )  where  the  upper  limit  is  the  value  of  z'  where 

is  reached.  (If  the  flood  plain  wall  shear  is  measured, 

could  also  be  computed  by  considering  all  the  forces 
acting  on  the  flood  plain  flow  for  unit  length.)  The  results 
for  all  the  experiments  are  given  in  Table  1.  Neglecting  the 
experiments  with  small  values  of  d  (since  in  such  cases,  the 
error  in  TQ  measurements  will  be  large),  one  general  comment 
that  could  be  made  is  that  Tq  is  of  the  same  order  as  the 
maximum  shear  stress  in  the  flood  plain  which  occurs  near 
the  junction  plane. 
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Using  the  results  presented  earlier,  it  is  possible  to 
present  an  analysis  of  the  velocity  profiles  in  the  main 
channel.  With  reference  to  Figure  3.1(d),  for  any  value  of 
z,  let  curve  1  represent  the  profile  of  u  with  y  that  would 
exist  when  the  side  wall  of  the  main  channel  is  extended  to 
a  depth  of  D.  If  now  the  extension  of  the  side  wall  is 
removed,  the  actual  velocity  profile  has  the  shape  of  curve 
3.  Curve  3  predicts  smaller  velocities  in  the  neighbourhood 
of  the  bed  because  of  the  reduction  of  T  .  Further,  a  y 
approaches  h  and  for  y>h  ,  due  to  the  direct  momentum  loss 
to  the  flood  plain,  the  velocity  profile  shows  a  dip  r 
defect  from  curve  2,  which  is  different  from  curve  1  because 
of  only  the  reduction  in  T*.  For  the  flood  plain,  in  the 
region  affected  by  the  interaction,  if  curve  4  in  Figure 
3.1(d)  is  the  undisturbed  velocity  distribution,  curve  5 
represents  a  typical  velocity  profile  which  includes  an 
increase  due  to  the  momentum  transport  to  the  flood  plain. 

It  was  found  that  for  a  number  of  the  experiments,  the 
shaded  area  between  curves  2  and  3  integrated  for  the  whole 
width  of  the  main  channel  was  approximately  equal  to  the 
shaded  area  between  curves  4  and  5  integrated  for  the  whole 
width  of  the  flood  plain.  That  means  the  area  between  curves 
1  and  2  integrated  to  cover  the  full  width  of  the  main 
channel  will  give  the  reduction  in  carrying  capacity  of  the 
compound  cross-section  because  of  the  flow  interaction.  It 
is  possible  to  express  this  observation  analytically. 
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3  .  4  Effect  o  f  Roughness 


Two  experiments,  experiment  9  and  10,  were  performed 
with  the  flood  plain  roughened  with  a  dense  hemispherical 
roughness  of  diameter  of  0.5  inch  (12.7  mm).  The  equivalent 
sand  roughness  of  this  roughness  has  earlier  been  found 
equal  to  0.028  feet  (8.5  mm)  (Ho llingshead ,  1972).  The 

results  of  the  velocity  measurements  showed  that  the 
momentum  transport  exists  as  in  the  case  of  a  smooth  bed. 

The  velocity  profiles  in  the  main  channel  (for  y>h)  and  in 
the  flood  plain  were  similar  when  considered  relative  to  the 
undisturbed  flood  plain  velocity  and  agreed  with  the  results 
for  the  smooth  bed. 

Considering  the  velocity  scale,  Figure  3.6(a)  shows 
that  the  variation  of  (um  -  u00)/(uqo)  with  y'/d  is  not  at 

all  affected  by  the  roughness.  The  length  scale  b  /d  is 

m 

increased  by  the  roughness  (see  Figure  3.7(a)).  The  velocity 

scale  (um  -  u00)/(u00)  is  affected  to  a  certain  extent 
(Figure  3.9(a))  and  the  length  scale  b^/d  is  considerably 
reduced  by  the  roughness  (Figure  3.9(c)).  It  is  necessary  to 
perform  more  experiments  with  a  range  of  roughness  heights 
before  the  effects  of  the  flood  plain  roughness  on  the 
momentum  transport  could  be  properly  assessed. 


3  .  5  Conclusions 

Based  on  the  experimental  results  presented  in  this 
chapter,  the  following  conclusions  could  be  made.  The 
transport  of  the  longitudinal  momentum  from  the  main  channel 
to  the  flood  plain  as  well  as  the  existence  of  apprecieble 
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shear  stresses  in  the  junction  plane  between  the  main 

channel  and  flood  plain  have  been  established.  The  lateral 

distribution  of  the  longitudinal  velocity,  in  the  upper  part 

of  interaction  region  between  the  main  channel  and  the  flood 

plain  if  viewed  with  respect  to  the  proper  velocity,  has 

been  found  to  be  separately  similar  in  the  main  channel  and 

flood  plain.  In  the  main  channel  interaction  region,  the 

velocity  scale  is  the  difference  between  the  maximum 

velocity  u^  in  the  main  channel  and  that  in  the  junction 

plane  (u^)  and  the  base  velocity  (for  similarity 

considerations)  is  u' .  The  length  scale  b  is  the  lateral 

m  m 

distance  from  the  centerline  where  (u  -  u')=0.25(u  -  u '  )  • 

m  m  m 

The  dimensionless  velocity  distribution  is  described  by 

equation  3.1.  The  dimensionless  velocity  (u  -  u  ) /u 

m  oo  oo 

varies  with  D/d  and  y'/d  (see  Figure  3.6).  Futher, 

(um  ”  uoo)  /uoo  var^es  mainly  with  D/d  and  to  a  small  extent 
with  y'/d  for  y'/d>0.35.  The  dimensionless  length  scale 
b  /d,  varies  mainly  with  D/d  and  futher  this  variation  is 
linear  and  is  described  by  equation  3.5. 

For  the  interaction  region  in  the  flood  plain,  the 
velocity  scale  is  (u^  “  u^)  and  the  length  scale  is  b^,  the 
distance  from  the  junction  plane  to  the  point  where 

“  Uoo  )  =0 . 5  (u^  ~  Uoo^*  The  dimensionless  length  scale  b^/d 
has  been  found  to  vary  mainly  with  D/d  and  this  variation 
is  described  by  equation  3.6.  In  the  main  channel  below  the 
level  of  the  flood  plain,  the  velocity  profiles  have  been 
found  to  be  approximately  similar  with  uq  the  velocity  in 
the  main  channel  at  y=h  as  the  velocity  scale.  The  value  of 
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uq  could  be  obtained  from  the  similarity  profiles  presented 
earlier. 

The  increase  in  the  shear  stress  in  the  flood  plain  in 

terms  of  x_  as  well  as  the  extent  of  the  affected  region 
0  oo 

in  terms  of  d  appear  to  be  mainly  function  of  D/d.  Most  of 
the  results  in  this  chapter  are  concerned  with  smooth 
boundaries  and  hence  more  experiments  should  be  performed 
with  different  kinds  and  heigths  of  roughness  elements 
before  the  results  presented  in  this  chapter  could  be 
adopted  for  practical  use. 


4.  Interaction  Between  Straight  Main  Channel  and  Straight 
Flood  Plain  Flows-Fur ther  S tudies-Asymme t rical  Channel 
In  Chapter  III  the  results  of  the  interaction  between  the 
main  channel  and  the  flood  plain  flows,  for  an  idealized 
case  of  a  straight  channel  in  a  straight  valley,  with  the 
main  channel  being  narrow  so  that  the  effects  of  interaction 
dominated  the  entire  main  channel,  was  studied.  In  this 
chapter  the  main  channel  is  wide  enough  to  have  at  least  a 
section  at  the  center  in  which  the  effects  of  lateral 
transfer  of  momentum  are  negligible.  Hence,  the  results 
presented  in  this  chapter  have  greater  potential  for 
applicability  to  practical  problems. 

In  this  chapter  mainly  the  effect  of  the  width  of  the 
channels  on  the  structure  of  the  flow  is  studied.  Although, 
three  types  of  asymmetrical  channels  (Type  A,B  and  C)  were 
studied,  greater  number  of  runs  were  made  with  type  C 
because  of  its  wide  ma in  channel  and  flood  plain  and  this 
chapter  presents  the  results  on  only  the  type  C  channel. 


4  .  1  Exp  eriments 

The  measurements  were  made  in  the  middle  cross-section 
of  the  channel  at  a  distance  of  approximately  30  feet  (9.15 
m)  from  the  entrance  after  making  sure  the  flow  was  fully 
developed  and  that  the  depth  of  flow  was  constant.  It  was 
believed  that  uniform  flow  was  established.  This  was  accom¬ 
plished  by  manipulating  the  tailgate  and  observing  that  the 
differential  pressure  reading  between  each  consecutive  pair 
of  the  seven  static  probes  was  the  same  up  to  one  millivolt 


65 


■  ’ 

F  ra[fl|.  'B2E!  *  /S  Kt  5*1  •mm  •  j 


66 


reading  on  the  modulator  of  the  transducer.  Also,  care  was 
taken  to  make  sure  the  water  'surface  slope  and  the  bed  slope 
were  parallel.  Five  experiments  were  conducted  in  the  type  C 
channel,  with  different  values  for  d  or  D.  Table  2  shows  the 
pertinent  details  of  the  dimensions  and  the  experiments. 

The  velocity  and  shear  stress  measurements  were  made  in 
the  whole  cross-section  of  the  compound  channel.  The 
boundary  surface  was  generally  smooth.  The  interaction 
regions  were  studied  first.  It  was  found  where  the 
interaction  would  take  place  and  this  area  was  studied  in 
detail.  The  vertical  velocity  profiles  in  the  mixing  zone 
were  very  close  together.  In  addition,  one  special 
experiment  was  conducted  with  the  multi-level  flood  plain 
and  because  of  its  special  features  the  results  of  this 
experiment  will  be  discussed  in  Appendix  B. 


4  .  2  Experimental  F  e  s  ul t  s 

The  velocity  profile  of  a  typical  experiment  are  shown 
in  Figure  A. 2.  Since  the  Froude  number  of  all  the 
experiments  was  around  0.35,  it  was  decided  to  present  the 
velocity  profiles  of  one  experiment  (number  3  with 
D/d=4.76).  It  was  found  that  the  velocity  will  continuously 
increase  as  one  moves  from  the  flood  plain  towards  the  main 
channel  and  reaches  a  point  where  the  disturbed  velocity 
approaches  the  local  undisturbed  velocity.  Figure  4.2  shows 
clearly  this  location  of  disturbed  velocity  approaching  the 
undisturbed  velocity.  In  the  main  channel  the  velocity 
increases  with  y,  up  to  some  level,  and  then  decreases 
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TABLE  2 


Significant  Details  of  the  Experiments 


B  =  2.35  feet,  b  =  1.67  feet 


Exp  t  . 

hCft) 

D  (  f  t  ) 

d(ft) 

0 

S 

D/d 

B/D 

b  /d 

No  . 

(cf  s ) 

x  1  0  3 

1 

0.323 

0.373 

0.050 

0.96 

0.450 

7.46 

6.31 

33.3 

3 

0.312 

0.395 

0.083 

1 . 06 

0.480 

4.76 

5 .95 

20.1 

4 

0.322 

0.4  76 

0.154 

1  .56 

0 .450 

3.09 

4.^4 

10.8 

5 

0.321 

0.5  96 

0.275 

1.98 

0.364 

2.17 

3 .94 

6.07 

6 

0.317 

0.353 

0.036 

0 . 90 

0.408 

9.81 

6.66 

4  6.4 

7* 

0 . 3  IP 

0.420 

0.102 

1.26 

0.724 

4.12 

5.61 

16.3 

8  * 

0.320 

0.500 

0. 180 

1.97 

0.600 

2 .73 

4.71 

9.26 

Dye  Test 


* 
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(a)  Sectional  View 


Figure  4.1  Definition  Sketches 
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Figure  4.2  Velocity  Profile:  Expt-3 
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continuously  to  eventually  become  approximately  constant 

near  the  water  surface.  This  behavior  was  also  observed  in 

the  previous  chapter  (Symmetrical  Channel).  In  the  central 

region  of  the  main  channel,  u  increases  continuously  with  y 

up  to  (almost)  the  water  surface  whereas  in  the  region  near 

the  flood  plain  interface,  a  dip  in  the  velocity  profiles 

due  to  lateral  momentum  transfer  could  be  seen.  At  least  in 

principle,  this  dip  is  similar  to  the  dip  in  velocity 

profiles  observed  in  narrow  channels  due  to  the  late1*-’..! 

momentum  transfer  to  the  side  wall  regions.  In  the  flood 

plain  region,  the  velocity  increases  continuously  with  y' 

until  the  water  surface  is  reached  except  in  the  narrow 

region  near  the  side  wall  of  the  flood  plain. 

Although,  velocity  profiles  in  the  lateral  direction 

(u(z))  in  the  region  where  y>h  for  a  few  values  of  y'  (the 

normal  distance  from  the  bed  of  the  flood  plain)  were 

prepared  for  all  experiments,  only  three  of  those  plots  are 

shown  in  Figure  4.3,  for  three  typical  experiments.  In 

Figure  4.3  for  the  three  experiments,  it  is  seen  that  u 

decreases  continuously  from  a  maximum  value  of  u  at 

m 

approximately  the  centerline  of  the  main  channel  to  u^  at 

the  interface  between  the  main  channel  and  the  flood  plain 

and  u  further  decreases  with  z'  to  reach  u  for  some  value 

oo 

of  z'.  The  value  of  u  remains  constant  at  u  for  a  portion 

oo  r 

of  the  flood  plain  until  the  section  where  influence  of  the 
wall  is  felt.  In  the  main  channel  for  z<0.0,  the  u  velocity 
decreases  due  to  the  influence  of  the  side  wall.  In  the 
immediate  vicinity  of  the  junction  plane,  in  most  of  the 
experiments,  a  sharp  increase  in  velocity  was  noticed.  This 
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Figure  4.3  (cont.)  Lateral  Velocity  Profiles:  Expt-3 
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sudden  increase  has  also  been  reported  by  other  observers 

such  as  Goncharov  (1964). 

Profiles  of  the  boundary  shear  stress  to  for  the 

experiments  are  shown  in  Figure  4.4.  Considering  the  bed 

shear  stress,  it  decreases  from  a  maximum  value  of  t  „  in 

om 

the  central  region  of  the  main  channel  to  some  value  as  the 
junction  plane  with  the  flood  plain  is  reached.  In  the  flood 
plain,  tq  decreases  from  a  large  value  at  z'  =  0.0  until  tqQ0 
is  reached.  The  presence  of  the  side  walls  causes  a  further 
decrease  in  the  bed  shear  stress. 


4 . 3  Analy s is  o  f  F  e s  u 1 1  s 

As  indicated  by  an  earlier  study  with  a  narrow  main 

channel,  u  velocity  profiles  in  the  main  channel,  in  the 

region  for  y>h,  in  the  lateral  direction,  (i.e.  u(z))  were 

plotted  with  (u  -  u')/(u  -  u')  versus  n=z/b  wherein  b  is 

m  m  m  1  m  m 

the  length  scale  chosen  in  such  a  way  that  at  j]  ~  1  •  0 , 

(u  -  u')/(u  -  u')=0.25.  The  results,  shown  in  Figure  4.5, 

m  m  m 

indicate  that,  in  general,  for  y'/d=0.0  to  1.0,  the  velocity 
profiles  are  similar  and  that  they  are  well  described  by  the 
equation: 


f  (n) 


U  -  U  '  o 

- —  =  1  -  0 . 7  5  n  2 

u  -  u ' 
m  m 


(4.1) 


Considering  the  velocity  scales,  if  the  main 
flood  plain  were  so  wide  that  sections  where 
profiles  (u(y)  or  u(y'))  are  not  affected  to 
(noticeable)  extent  by  the  momentum  transfer 


channel  and 
the  velocity 
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to  the  flood 
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Figure  4.4  Bed  Shear  Stress:  Expts  1  &  3 
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Figure  4.4  (cont.)  Bed  Shear  Stress:  Expts  4  &  5 
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Figure  4.4  (cont.)  Bed  Shear  Stress:  Expts  6,  7  &  8 
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8  0 


plain  and  assuming  that  in  those  unaffected  regions,  the 
vertical  velocity  profiles  (u(y)  or  u(y'))  are  described  by 
the  logarithmic  law,  then: 


u  -  u 
m  oo 


u 


co 


C  .  D 

*m 

C*f  d 


[  1  + 


log 


C 

log  — 


D 

y'/c*t  d 

y '/ C*f  gds 


1  (4.2) 


wherein  C  is  a  constant  equal  to  9.06  and  C.  =  T  /  y  D  S  and 

n  *m  om'  '  o 

C+f=Tn  /YdS  .  The  factors  C+  and  could  be  obtained 

*f  Ooo  o  *m  *f 

from  Figure  A- 2  in  Appendix  A.  For  the  present  experiments, 

values  of  ( u m  -  Uoo^^u00)  calculated  from  equations  4.2  are 
compared  with  measured  values.  Figure  4.6  shows  that  the 
calculated  values  agree  reasonably  well  with  the  measured 
values. 

To  predict  uffi,  (um  -  uoo)/(um  -  u^)  was  plotted 
against  y'/d  for  the  seven  experiments.  It  appears  from 

Figure  4.7  that  (um  “  Uco^  ^Um  ~  Uoo^  varies  mainly  with  y'/d 
and  increases  from  about  0.2  for  y'/d=0.0  to  about  0.30 
for  y'/d=1.0.  An  average  value  of  0.25  could  be  suggested 
for  practical  purposes. 

The  data  for  the  length  scale  showed  that  for  any 

experiment,  the  variation  of  b^  with  y'  was  small  hence 

average  values  for  b  /D  for  the  whole  range  of  y'/d  were 

m 

obtained.  These  values  were  then  plotted  against  D/d  in 
Figure  4.8.  From  Figure  4.8  it  is  seen  that  b^/d  increases 
linearly  with  D/d  and  could  be  described  by  the  equation: 


S' 
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(meas) 
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Figure  4.6  Calculated  Velocity  vs.  Measured  Velocity 
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Figure  4.7  Variation  of  Normalized  Depth  vs.  Normalized  Velocity 


83 


Figure  4.8  Length  Scales:  (bf,  bm,  bt,  br) 


84 


b  n 

— ^  =  3 . 7  8  (  —  -  1  ) 

d  d 


(4.3) 


The  velocity  profiles  in  the  flood  plain  are  tested  for 

similarity  by  plotting  (u  -  u  )/(u'  -  u  )  versus  z'/b£ 

UJ  m  oo  f 

where  bf=z'  where  (u  -  u00)/(u^  -  Uqo)=0.50.  Figure  4.9 
shows  that,  in  general,  these  profiles  are  similar  and  could 
be  described  by  the  exponential  equation: 


f  (tT) 


u-u  -0.693  n  "  ^ 

00  1 

7 -  =  e 

u  -  u 
m  oo 


(4.4) 


where  r)'=z'/b^.  The  undisturbed  velocity  in  the  flood  plain 

could  be  predicted  using  the  logarithmic  law  and  then  u  ' 

m 

could  be  calculated  using  the  results  presented  earlier. 
Figure  4.8  shows  that  the  dimensionless  length  scale  b^/d 
increases  linearly  with  D/d  and  could  be  described  by  the 
equation: 


— T-  =  0 . 6  4  (  —  -  1  )  (4.5) 

d  d 

The  width  of  the  momentum  diffusion  region  could  be 
seen  to  start  from  approximately  the  section  in  the  main 
channel  where  u  begins  to  decrease  with  z  to  the  section  in 
the  flood  plain  where  u  reaches  .  Using  equations  4.1, 

4.3,  4.4  and  4.5  the  total  width  of  the  length  scale  b  ^ ,  was 

calculated.  The  results  plotted  in  Figure  4.8  show  that  b^/d 
increases  linearly  with  D/d  and  is  described  by  the 


equation : 


, 
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Figure  4.9  Velocity  Profile  Similarity:  Flood  Plain  (Expts  1,  3,  4,  5  &  6) 
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(4.6) 


The  dotted  line  in  Figure  4.8  shows  the  width  of  the 
diffusion  region  from  an  earlier  study  (Chapter  III)  where 
the  main  channel  was  so  narrow  that  it  did  not  allow  the 
full  growth  of  the  diffusion  region. 

The  experimental  results  shown  in  Figure  4.4  indicated 
that  the  bed  shear  stress  in  the  flood  plain,  was  increased 
due  to  the  interaction.  Figure  4.10  shows  that  the  profiles 
of  the  excess  bed  shear  stress  in  the  flood  plain  defined  as 

AV(I°  '  I0oo)  are  silcllar  when  AT0/(Ton  -  T0co)  ls  Plotted 

against  z  ' /b  T  where  b  =z '  where  /  (t'  _  -  x  «  )  =0 . 5 .  The 

1  T  o  om  Ooo 

similarity  curve  is  well  represented  by  the  equation: 


f  (  n') 
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Ato 


T  '  -  T 

om  Ooo 


,2 

-0.693nT 

e 


(4.7) 


The  behaviour  of  the  scale  for  the  excess  bed  shear  stress 
(T^m  -  T )  /  (T0oo )  and  the  length  scale  bT  /d  is  shown  in 
Figure  4.10.  Eoth  these  scales,  when  expressed  in  a  dimen¬ 
sionless  manner,  increase  linearly  with  D/d  and  could  be 
represented  by  the  equations: 


T'  _  T 

om  Ooo 
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0  oo 
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(4.8) 
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(4.9) 


* 


87 


Figure  4.10  Bed  Shear  Stress  Similarity:  Flood  Plain  Channel  (Expts  1  to  8) 
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Considering  the  bed  shear  stress  in  the  main  channel, 

it  is  known  that  it  is  reduced  due  to  the  transfer  of 

moment  urn  for  the  flood  plain.  The  results  shown  in  Figure 

4.11  support  this  idea.  If  Tom  is  the  bed  shear  stress  at 

approximately  the  middle  of  the  main  channel,  wh i c h  is 

apparently  not  affected  by  the  lateral  momentum  transfer, 

and  Atq  is  the  reduction  in  shear  stress  at  any  z  due  to  the 

momentum  transfer.  Figure  4.11  shows  that  At  / i"  reases 

with  =z/B/2  for  each  value  of  D/d.  If  (At  /T  )  is  the 

o  on  n 

maximum  value  of  the  reduction  as  obtained  from  Figure  4.11, 
it  is  seen  that  (ATo/Tom)^  i s  mainly  a  function  of  D/d.  For 
a  D/d  =  10.0,  (At  / T  )  is  about  0.20. 

o  om'm 

It  was  questioned  as  to  whether  the  velocity 
distribution  in  the  region  affected  by  the  lateral  momentum 
transfer  is  still  logarithmic.  From  Figure  4.12  for 
experiment  1  and  other  plots  (not  reproduced  herein)  it  was 
found  that  in  the  portions  of  the  main  channel  and  flood 
plain  that  were  not  affected  by  the  momentum  transfer,  the 
velocity  distribution  follows  the  logarithmic  law: 


u 

m 


y  u* 

5.75  1 o  g (  -  )  +  5.5 


(4.10) 


except  in  the  narrow  region  close  to  the  side  walls.  It  was 
also  found  that  in  the  interaction  region,  equation  4.10 
describes  the  data  well  if  the  local  value  of  u^  is  used  in 
the  equation.  This  observation  suggests  a  possible 
analytical  attack  on  the  interaction  region  combining  the 
law  of  the  wall  and  shear  layer  concepts.  In  the  main 
channel,  in  the  interaction  region,  the  momentum  transfer  to 
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Figure  4.1 1  Variation  of  Bed  Shear  Stress  with  D/d  &  X 
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Figure  4.12  Velocity  Distribution  (u(y')):  Log  Law 
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Expt  -  1  (c) 
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Figure  4.12  (cont.)  Velocity  Distribution  (u(y)):  Log  Law 
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Figure  4.12  (cont.)  Velocity  Distribution  (u(y) ) :  Log  Law 
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the  flood  plain  causes  a  significant  dip  (more  noticeable 
than  that  in  channels  without  a  flood  plain)  in  the  velocity- 
profile  in  the  upper  part  of  the  flow.  If  y*  ±s  the  location 
of  the  start  of  the  dip  in  the  velocity  profile  (see  Figure 
4.12)  in  Figure  4.13  y*/D  is  plotted  against  A=z/B/2  for  a 
number  of  experiments.  One  can  at  most,  indicate  some  trends 
in  Figure  4.13.  It  appears  that  for  any  given  D/d,  y*/D 

JL 

decreases  as  X  increases  and  for  X  close  to  unity,  y  /d 
attains  values  in  the  range  of  0.45  to  0.7. 


4 . 4  Kinema  tic  F  d  dy  Viscosity 

It  was  shown  that  the  vertical  distribution  of  the  ve¬ 
locity  even  in  the  region  affected  by  the  lateral  momentum 
transfer  is  logarithmic.  Average  kinematic  eddy  viscosity  in 
y  direction  £  ,  could  then  be  shown  to  be  equal  to 

y 

0.067  u  ^  d . 

Dye  injection  was  adopted  to  find  out  the  magnitude  and 

variation  of  the  lateral  kinematic  eddy  viscosity  £  a"  o 

z 

see  whether  it  was  affected  by  the  lateral  momentum  transfer 
in  the  interaction  region.  It  can  be  shown  the  that  rate  of 
growth  of  a  passive  jet  in  a  media  with  the  same  density  as 
the  fluid  injected  is  a  function  of  £  ,  x  and  u.  The  final 
working  equation  is  given  (lecture  notes,  Pajaratnam  1978) 
as : 

/£  X 

— —  (4.11) 

where  b  is  the  half  width  of  the  plume  at  distance  x  from 
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Figure  4.13  Location  of  Dip  in  Velocity  Profiles:  Main  Channel 
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the  nozzle  and  u  is  the  velocity  at  the  point  where  the 
nozzle  is  located. 

Lateral  distribution  of  the  velocity  was  defined  by 
equations  4.1  and  4.4  in  the  main  channel  and  flood  plain 
channel  respectively.  Therefore,  could  be  calculated  as 

f  o 1 lows : 


T 


o  z  x 


(4.12) 


Equation  4.12  could  be  evaluated,  since  the  values  of  T 

o  zx 

could  be  obtained  simply  by  considering  the  equilibriun 
condition  for  a  column  of  water  (see  sketches  in  Figure 
4.14). 

Two  series  of  dye  test  were  conducted  in  the 

asymmetrical  channel  Type  C.  The  dye  had  a  density  equal  to 

the  water  and  was  injected  through  a  one  milimeter  OD  tube. 

The  dye  was  injected  parallel  to  the  x  direction  at  half 

depth  on  the  flood  plain.  A  constant  head-tank  provided  the 

necessary  head  in  such  a  way  that  the  dye  could  flow  with 

the  same  velocity  as  the  surrounding  flow.  The  velocity  u 

and  bed  shear  stress  tq  were  measured  for  each  profile.  Each 

time  the  dye  was  injected,  a  few  photos  were  taken  from  four 

feet  (1.25  m)  above  the  water  surface.  The  photos  provide  a 

measure  of  the  rate  of  expansion  of  the  jet. 

From  equation  4.11  the  value  of  lateral  eddy  viscosity 

£  was  calculated.  The  values  of  £  ,  and  £  at  the  level  of 
z  z  z 

y'=d/2=0.05  feet  and  y  '  =  0 . 0  0  5  feet  (1.5  and  0.15  cm)  for 
experiment  number  7  and  y'=d/2=0.09  feet  (2.75  cm)  for 
experiment  number  8  were  obtained  from  the  dye  test.  Figure 
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Figure  4.15  Cross-Sectional  Isovel  Map:  Expt-4 
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4.14  shows  the  lateral  distribution  of  £  and  £  ,  for 

z  z 

experiments  with  one  large  and  one  small  value  of  D/d 

(experiment  number  1,  D/d=7.46  and  experiment  number  5, 

D/d=2.17).  Experiment  number  7  is  also  included. 

Variation  of  the  transverse-mixing  coefficient  £  /d  u  .  , 

over  the  width  of  the  channel  is  shown  in  Figure  4.14. 

Fischer  (1973)  showed  that  the  transverse-mixing  coefficient 

is  a  function  of  the  aspect  ratio.  For  a  range  of  aspect 

ratios  from  3  to  60,  the  transverse -mixing  coefficient 

varies  from  0.09  to  0.20.  For  each  set  of  experiments,  the 

aspect  ratio  is  constant.  Therefore,  the  variation  of  the 

transverse-  mixing  coefficient  over  the  width  of  the  channel 

in  Figure  4.14  shows  the  effect  of  the  interaction  between 

the  main  flow  and  the  flood  plain  flow  on  £  .  It  is  shown 

z 

that  in  the  mixing  region,  £  or  £  ,  is  increased.  The 

z  z 

transverse-  mixing  coefficient  in  the  interaction  region 
decreases  rapidly  from  the  junction  plane  toward  the  flood 
plain  bank  and  it  decreases  slowly  toward  the  center  of  the 
main  channel.  Depending  on  the  D/d,  the  transverse -mixing 
coefficient  in  this  mixing  region  will  assume  appreciably 
larger  values,  more  than  the  undisturbed  value  for  a  channel 
with  the  same  aspect  ratio  of  the  flood  plain  or  main 
ch  annel . 

A  cross-sectional  isovel  map  of  experiment  number  4  is 
shown  in  Figure  4.15.  From  this  isovel  map,  one  can 
visualize  the  zone  of  interaction,  the  direction  of  the 
momentum  transfer  and  the  change  in  the  magnitude  of  the 
veloci ty . 
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4 . 5  Con  clu  s ions 

The  results  presented  in  this  chapter  demonstrate  that 

the  structure  of  the  flow  in  the  main  channel  with  a  flood 

plain  could  be  analysed  using  the  concept  of  lateral 

momentum  transfer  and  by  treating  the  interaction  region  as 

a  turbulent  shear  layer.  The  interaction  between  the  main 

channel  and  the  flood  plain  flow  produces  increased 

velocities  and  bed  shear  stresses  in  the  flood  plain  and 

causes  reductions  of  velocities  and  bed  shear  stress  in 

the  main  channel.  These  changes  have  been  analysed  using 

similarity  principles.  It  has  been  shown  that  in  the  upper 

part  of  the  flow  (y>h),  in  the  main  channel  and  in  the  flood 

plain  area  of  the  interaction  zone,  the  lateral  velocity 

profiles  are  similar  and  can  be  described  by  equations  4.1 

and  4.^  The  related  length  scales,  b  and  b  in  terms  of  d, 

f  m 

the  depth  of  flow  in  the  flood  plain  increase  linearly  with 
D/d,  which  are  defined  by  equations  4.3  and  4.5.  Further 
findings  reveal  that  the  velocity  profiles  in  the  vertical 
direction  are  described  by  the  logarithmic  velocity 
distribution  law  (equation  4.10)  if  the  local  bed  shear 
stress  is  used  as  the  velocity  scale,  except  in  the  upper 
part  of  the  flow  in  the  main  channel.  Equations  have  also 
been  developed  for  predicting  the  revelant  velocity  scales 
(equation  4.2  and  Figure  4.6  ). 

From  the  variation  of  the  maximum  bed  shear  stress  with 
aspect  ratio  shown  in  Figure  A-l,  and  the  lateral 
distribution  of  the  bed  shear  stress  in  the  undisturbed 
zones  (main  channel  and  flood  plain)  shown  in  Figure  A-2 » 
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the  variation  of  tq  in  the  interaction  region  in  the  main 
channel  could  be  predicted  (Figure  4.11).  Equation  4.8 
provides  the  magnitude  of  T'Qm  at  the  junction  plane.  The 
lateral  distribution  of  the  excess  bed  shear  stress  in  the 
flood  plain  is  presented  by  the  similarity  curve,  equation 
4.7.  The  length  scale  for  the  bed  shear  stress  in  the  flood 
plain  varies  linearly  with  D/d,  equation  4.9.  Proper 
manipulation  of  equations  4.7  and  4.9  provides  tq  in  the 
affected  region  of  the  flood  plain.  This  suggests 
possibilities  of  predicting  the  magnitude  of  the  turbulent 
eddy  viscosity  for  the  interaction  region  analytically  and 
this  procedure  is  given  in  Figure  4.14.  From  this  figure,  it 
could  be  seen  that  the  eddy  viscosity  in  the  junction  plane 
is  increased  sharply. 
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5.  Meandering  Main  Channel  with  Straight  Flood  Plains 
Although  several  studies  have  been  made  on  pipe  bends  and 
open  channel  bends,  only  a  few  investigators  have  studied 
meandering  channels  with  flood  plains.  Rajaratnam  & 
Muralidhar  (Personal  Communication)  in  1967  collected  a  list 
of  over  100  references  related  to  the  subject  of  open 
channel  curves  and  pipe  bends. 

The  fact  is,  because  of  the  presence  of  centripetal 
acceleration,  the  water  surface  is  superelevated;  spiral 
motion  is  set  up,  and  the  velocity  and  boundary  shear  stress 
distributions  will  be  modified  accordingly.  The  flow 
characteristics  will  change  from  section  to  section.  The 
effect  of  the  bend  extends  in  all  directions.  The  bend  is  an 
obstacle  to  the  flow,  causing  additional  resistance  besides 
the  boundary  resistance.  As  a  result,  there  is  extra  energy 
loss  and  a  backwater  curve  is  created.  In  tight  bends, 
separation  may  occur.  In  natural  rivers,  the  flow  is  further 
complicated  by  irregular  channel  geometry  and  movable  bed 
material.  The  presence  of  the  flood  plain  brings  additional 
complications  to  the  study  of  momentum  exchange  between  the 
flow  in  the  meandering  river  and  the  flood  plain  flow.  A 
brief  literature  review  on  curved  was  presented  in  Section 
1.3.  channels • 

5 . 1  F low  i n  Curved  Channels 

A  good  review  of  the  nature  of  the  flow  in  the  curved 
channel  has  been  written  by  Rozovskii  (1961).  He  accounts 
for  the  existing  forces  acting  on  the  flow.  Because  of  the 
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curvature  of  the  channel,  the  streamlines  are  curved  in  plan 
and  this  results  in  transverse  inclination  of  the  free  water 
surface.  The  extent  of  this  superelevation  can  be  evaluated 
approximately  by  considering  the  equilibrium  conditions  for 
a  column  of  water  on  the  bed.  Actually,  the  three  main 
forces,  namely  the  centrifugal  force  F  ,  differential 
pressure  on  the  sides  of  the  prism  AP  ,  and  frictional  force 
Ff  exerted  on  the  column  of  fluid  must  be  in  equilibrium  and 
using  this  observation,  one  can  derive  an  expression  for 
radial  slope  S^.  If  frictional  force  is  ignored  one  could 
show  that: 
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(5.1) 


An  expression  for  transverse  circulation  can  be  derived 
by  considering  equilibrium  for  a  volume  element  assuming  the 
radial  component  can  be  Ignored. 


V 

Sr  =  -  (5.2) 

r  g  r 

Since  equation  5.2  is  independent  of  z,  this  equation  can 
only  be  satisfied  in  two  cases: 

1.  Straight  channel,  r=infinity  and, 

2.  constant  velocity,  V=constant.  Because  of  bed  friction, 
the  vertical  velocity  varies  from  a  maximum  at  the  top 
and  a  minimum  value  at  the  bed.  The  only  place  equation 
5.2  can  be  satisfied  is  where  the  velocity  is  equal  to 
the  average  velocity. 

a.  For  particles  moving  with  a  higher  velocity  than  the 
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average  velocity,  the  centrifugal  forces  will  exceed 
the  differential  pressure.  Therefore,  these 
particles  will  be  displaced  in  a  radial  direction 
away  from  the  center  of  curvature, 
b.  In  the  lower  part  of  the  flow,  where  the  velocities 
are  lower  than  the  average  velocity,  the  centrifugal 
force  is  smaller  than  the  AP  and  particles  will 
move  toward  the  center  of  curvature. 

Roth  transverse  slope  and  helicodial  circulation  will 
change  the  velocity  structure  of  the  stream  bend.  The 
velocity  components  in t ens i f y  the  exchange  o f  momentum  in 
all  directions  especially  in  the  longitudinal  and  radial 
direct  ions  in  turbulent  flows  wh  ich  leads  t  o  a_ 
redistribution  o f  vert ica 1  and  rad ial  ve lo c i t ie s  .  This 
greater  change  o  f  momentum  over  the  width  o  f  the  curved 
ch anne Is  in  meandering  channels  with  the  flood  plain  channel 
adds  an  additional  comp 1 i cat  ion  t  o  the  exchange  o  f  momentum 
between  the  main  channel  and  the  flood  plain. 

A  brief  discussion  and  analysis  of  the  properties  of 
streams  in  bends  and  a  few  available  models  are  given  in  the 
following  section. 


5 . 2  Theoretical  Con  s  ide  ra  t  ions  o  f  Flow  i  n  Bend 

If  a  suitable  cylindrical  coordinate  system  is  chosen, 
the  Reynolds  equation  of  motion  and  the  continuity  equation 
theoretically  will  give  adequate  base  to  solve  the  problem 
of  flow  in  a  bend.  It  is  very  unlikely  that  a  complete 
solution  to  these  equations  of  motion  can  be  found.  For  an 
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approximate  solution  to  this  problem,  a  simplified  form  of 
these  equations  usually  based  on  the  assumptions  of  fully 
developed  flow,  (see  Fozovskii  1961,  B.C.  Yen  1965,  de 
Vriend  1976,)  flat  bed,  and  or  constant  loss  of  energy  for 
an  equal  length  of  flow,  could  be  established.  Un f o r t u na t ly , 
even  in  this  simplified  form,  there  are  ten  unknowns  in 
those  four  nonlinear  differential  equations.  Therefore,  at 
least  six  more  relations  are  necessary.  Even  with  these 
conditions,  a  complete  solution  to  such  a  highly  nonlinear 
differential  equation  would  still  be  impossible.  Therefore, 
further  approximations  and  simplification  have  to  be 
introduced.  The  latest  simplification  brings  restrictions  on 
the  geometry  of  the  channel.  The  distribution  of  the 
tangential  velocity  components  u^  has  to  be  assumed  ahead  of 
time.  Restrictions  are  as  follows: 

1.  geometrical  restrictions,  which  state  the  applicability 

of  the  model  at  the  following  points: 

a.  the  channel  must  be  shallow  which  means  a  large 
aspect  ratio  (B/D>10.0); 

b.  width  of  the  channel  and  the  radius  to  be  of  the 
same  order  of  magnitude; 

2.  flow  restrictions: 

a.  at  the  core  of  the  flow,  all  vertical  profiles  will 
assume  a  predetermined  velocity  distribution;  most 
investigators  are  suggesting  the  logarithmic 
velocity  profile  for  their  models; 

b.  similarity  assumption  for  the  distribution  of  the 
lateral  velocity  components  is  valid; 

c.  the  shear  stress  effects  must  dominate  the  inertial 
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effects. 

F.  Engelund  (1975)  introduced  his  simple  model  by  providing 
a  relationship  based  on  the  above  assumptions.  His  approach 
is  as  f  o  Hows  : 

1.  simplified  equations  of  motion  are: 
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(5.3) 


(5.4) 


where  h  the  level  of  water  surface; 

2.  geometry  condition: 

S  r  =  S  R  (5.5) 

o 

3.  the  slip  velocity  method  was  chosen  for  the  distribution 
of  Uq  over  the  depth.  The  logarithmic  law  was  assumed 
near  the  bed  and  the  parabolic  distribution  was  assumed 
for  the  region  away  from  the  wall  region. 

4.  Similarity  assumptions  are  as  follows: 
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where  U  is  some  reference  velocity. 

After  solving  for  the  pair  of  ordinary  differential 
equations  and  imposing  boundary  conditions  to  obtain  the 
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values  for  the  integration  constants,  Engelund  obtained 
expressions  for  f^and  f^« 

More  elaborate  models  are  those  of  Rozovskii  (1961), 
R.C.  Yen  (1965)  and  de  Vriend  (1976).  Eesides  the  complexity 
of  their  models,  they  do  not  offer  a  better  solution  of  the 
flow  around  the  bends.  Even  with  the  use  of  these  models, 
prediction  of  the  flow  does  not  match  the  measurements  and 
in  some  sections  the  trend  is  in  the  opposite  direction.  In 
chapter  I,  Section  1.3.2  some  remarks  were  made  on  the 
subject  and  it  was  indicated  why  analytical  solutions, 
considering  present  knowledge  of  flow  in  curved  channels, 
fail  to  predict  the  behavior  of  the  flow.  Also,  it  was 
mentioned  that  most  of  the  assumptions  necessary  to  simplify 
the  equations  of  motion  are  not  necessarily  true. 


5 . 3  Experiment  and  Exp  erimental  Results  on  Flow  i n  Curved 
C  h annels 

The  details  of  the  construction  of  the  curved  channel 
were  presented  in  section  2.3.3.  Chapter  II  also  described 
the  equipment  used  in  connection  with  the  measurement  of  the 
slope,  discharge,  velocity  and  shear  stress. 

The  measurements  were  made  at  the  beginning  of  the 
second  loop  at  a  distance  of  approximately  30  feet  (9.15  m) 
from  the  entrance.  Eighteen  feet  (5. 49  m)  of  the  straight 
main  channel  was  connected  with  the  first  loop  through  a 
gentle  curve.  It  was  believed  that  developed  flow  was 
obtained  at  the  exit  of  the  first  loop.  The  cross-sections 
were  located  at  15  degree  increments  of  9  from  the  first 
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section,  located  at  the  beginning  of  the  second  loop.  The 
sections  in  the  bend  were  named  by  their  angle  in  degrees 
and  also  numbered  from  1  to  15  where  1  corresponds  to 
entrance  section  and  15  corresponds  to  exit  section  of  the 
second  meandering  loop.  All  sections  were  perpendicular  to 
the  longitudinal  axis  of  the  curved  channel.  At  every 
cross-section,  nine  vertical  profiles  were  selected  in  such 
a  way  that  they  were  one  inch  (2.54  cm)  apart  from  each 
other.  These  profiles  were  named  by  the  distance  from  the 
inside  wall  of  the  curved  channel  (i.e.  R2  means  two  inches 
from  the  wall  or  a  distance  of  27  inches  from  the  center  of 
the  curve). 

One  experiment  with  depth  of  flow  equal  to  0.10  feet 

(3.05  cm)  was  performed.  Magnitude  and  the  angle  of 

deviation  of  the  velocity  from  the  tangential  direction  were 

measured  and  the  velocity  vector  was  resolved  into  the 

tangential  and  radial  components.  The  information  obtained 

at  any  point  such  as  x,y,  z,  R,  9,  V,  u  and  v  were  entered 

v  r 

into  a  matrix  for  further  processing  of  the  data  (i.e. 
plotting  the  isovels  etc.). 

Figure  5.1  shows  tangential  velocity  profiles  for 
cross-section  numbers  1  to  15.  Figure  5.2  presents  radial 
velocity  profiles.  Since  the  angle  of  deviation  was  small, 
the  tangential  velocities  and  the  velocity  vector  were 
almost  equal.  Therefore,  it  was  decided  not  to  plot  vector 
velocity  profiles.  A  quick  examination  of  tangential 
velocity  profiles  revealed  that  none  of  these  profiles 
followed  logarithmic  distribution.  However,  the  aspect  ratio 
for  this  run  was  about  8.3  (B/D=8.3)  which  was  not  large 
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Figure  5.1  Distribution  of  Tangential  Velocity  in  Curved  Channel,  Cross-Sections  1  to  15 
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Figure  5.2  Distribution  of  Radial  Velocity  in  Curved  Channel,  Cross-Sections  1  to  15 
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enough  to  neglect  the  effect  of  the  side  walls  on  the  core 
of  the  flow  but  it  was  similar  to  one  of  the  B.C.  Yen's 
experiments,  where  he  found  that  the  velocity  profiles  in 
the  center  of  the  channel  could  be  described  by  the  log-law. 
The  photos  in  Figure  5.3  shows  the  helicoidal  motion  of  the 
flow  in  the  curved  channel.  In  Figure  5.4  isovel  maps  of  the 
cross-sections  are  shown.  The  isovel  contours  of  any  given 
cross-section  are  based  on  approximately  200  data  points. 
Surface  II  program  (Sampson,  1978)  was  employed  for 
contouring  the  isovel  maps  and  the  Cal  comp  plotter  was  used 
for  drawing  the  isovel  maps.  Basically,  these  maps  we„a 
prepared  for  comparison  with  isovel  maps  of  the  main  channel 
with  the  flood  plain.  Lateral  distribution  of  the  depth 
averaged  tangential  velocity  is  shown  in  Figure  5.5.  It  can 
be  seen  from  the  figure  that,  these  average  tangential 
velocity  profiles  are  not  the  same.  It  can  easily  be  noticed 
that  the  maximum  velocity  at  the  beginning  of  the  bend  is 
close  to  the  inner  wall  and  gradually  shifts  toward  the 
other  side  (outer  wall).  By  the  time  the  flow  reaches  the 
end  of  bend  and  enters  the  next  bend,  the  order  of  shifting 
velocity  is  changed  and  again  the  maximum  velocity  will 
occur  near  the  inner  wall. 

Regarding  the  distribution  of  the  velocity  and 

applicability  of  the  defect  law  along  the  depth  of  flow,  the 

vector  and  tangential  velocities  were  plotted  and  compared 

with  a  theoretical  curve.  It  was  interesting  to  see  that 

none  of  these  profiles  obeyed  the  logarithmic  law.  The 

dimensionless  form  of  u^/u-  vs .  y/D  is  shown  in  Figure 

0  ©ma  x 

5.7.  As  can  be  seen  for  the  velocity  profile  located  in 


' 


' 


no 


Curved  Channel 


Helicoidal  Notion  in  the  Curved  Channel 


Helicoidal  Notion  in  the  Curved  Channel 
Figure  5.3  Photograph  :  Helicoidal  Motion 
of  the  Flow  in  the  Meandering  Channel 
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Figure  5.6 
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Figure  5.7  Variation  of  Velocity  and  Depth  Along  Centerline 
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Fluctuation  Of  u  &  D  Along  (£ 

Figure  5.8  Non-dimensional  Centerline  Velocity  Profiles  in  Curved  Channel,  Cross-Sections  1  to  15 

Lateral  u  Distribution 


Figure  5.9 


Legend  Same 
As  Fig.  5.7 


Non  dimensional  Lateral  Distribution  of  Tangential  Velocity  in  Curved  Channel, 
Cross-Sections  1  to  15 
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cross-section  number  1  (0=0.0  degree),  the  velocity 

increases  from  zero  (on  the  bed)  to  a  maximum  at  a  relative 
depth  of  0.2  (y/D=0.2).  Then  the  velocity  is  constant  up  to 

the  water  surface.  However,  in  cross-section  number  5  (0=45 
degrees)  the  velocity  increases  with  y  and  the  maximum 
occurs  approximately  near  the  water  surface.  Not  only  does 
the  vertical  distribution  of  tangential  velocity  vary  with 
0,  but  the  magnitude  of  the  velocity  depends  on  the  angle  0 
i.e,  the  location  of  cross-section  in  the  bend.  The  average 
velocity  of  each  cross-section  u  was  computed  for  all  the 
profiles.  If  the  average  velocity  of  cross-sections  1,  F  and 
15  which  correspond  to  0=0.0  degrees  is  set  to  U,  then  the 
variation  of  u/  U  with  0  for  cross-sections  1  to  15  is  shown 
in  Figure  5.8.  As  was  discussed  in  chapter  I,  the  slope  of 
the  main  channel  will  change  according  to  S  Cos0. 
Consequently,  with  the  increase  of  0  the  bed  slope  will  be 
reduced  and  therefore,  the  velocity  which  is  proportional  to 
the  water  surface  slope  to  some  power  will  be  decreased. 
Subsequently,  the  continuity  equation  reveals  that,  any 
reduction  in  the  average  velocity  requires  a  higher  average 
depth  of  flow  in  order  that  an  equal  amount  of  flow  passes 
through  all  cross-sections.  In  the  same  figure,  the 
variation  of  D/D^  (where  D  and  are  the  average  depths  of 

the  flow  at  any  section  and  that  at  0=0.0  degrees)  with  9  is 
shown.  Apparently,  the  reduction  of  the  velocity  at  0=45 
degrees  is  about  14  percent  (U  at  45  degrees/U  at  0=0.0  is 
equal  to  0.86)  and  the  increase  of  the  depth  of  flow  for  the 
same  cross-section  is  about  11  percent.  The  discrepancy  of 
3%  is  possibly  due  to  the  error  in  the  experimental 
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me  asur  ement  s . 

Non-dimensional  distribution  of  the  depth  averaged 

tangential  velocity  (  u/  u  where  u/  u  is  the  maximum 

J  max  max 

value  of  u)  over  the  width  of  the  channel  (r  -r.)/(r  -  r.), 

i  o  1 

where  r.,  r  and  r  are  innermost  radius,  outermost  radius 
i  o 

and  radius  of  the  bend  respectively,  is  shown  in  Figure  5.9. 
It  can  be  seen  from  these  lateral  velocity  distributions 
that,  at  cross-section  number  1  the  maximum  velocity  occurs 
near  the  inner  wall  and  the  lower  velocities  are  located 
near  the  outer  wall.  As  the  flow  goes  around  the  bend,  (i.e. 
as  9  increases)  the  maximum  velocity  will  shift  toward  the 
outer  wall.  Maximum  shifting  occurs  at  9=45  degrees  which  is 
half  of  the  central  angle  (9C=90  degrees).  All  of  these 
distribution  profiles  can  be  superimposed  on  top  of  each 
other  by  rotating  their  ordinates  around  a  point  located  at 
20  percent  of  the  depth  and  95  percent  of  the  maximum 
velocity.  Once  again,  the  importance  of  9  can  be  seen.  As 
Choudhary  &  Narasimhan  (1977)  and  Sieber  &  Gotz  (1975) 
showed  in  their  experimental  results,  the  flow  is  indeed 
dep  endent  on  9 . 

Generally,  radial  velocities  are  small  compared  with 
tangential  velocity.  The  maximum  radial  velocity  occurs  near 
the  bed  and  its  magnitude  is  about  10  percent  of  the  maximum 
tangential  velocity.  An  attempt  was  made  to  find  out  if  any 
of  the  similarity  hypotheses  suggested  by  previous  authors 
can  be  adopted  for  the  distribution  of  the  radial  velocity 
component.  It  seems  that  none  of  those  similarity  hypotheses 
hold  perfectly,  but  generally  the  prediction  of  radial 
velocity  and  the  measurements  have  the  same  trend. 
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The  variation  of  the  bed  shear  stress  in  the  meandering 
channel  can  be  seen  in  Figure  5.6.  The  maximum  bed  shear 
stress  is  located  near  the  inner  wall,  but  it  is  not 
symmetrical  with  respect  to  the  crest  of  the  curved  channel 
(cross-section  number  8).  The  maximum  bed  shear  stress  is 
located  close  to  the  inner  wall  and  it  can  be  as  large  as 
2.15  times  the  centerline  bed  shear  stress  of  the 
corresponding  straight  channel.  The  minimum  bed  shear  stress 
usually  occurs  near  the  outer  wall  and  is  only  60  percent  of 
the  centerline  bed  shear  stress  of  the  equivalent  straight 
channel . 


5 . 4  Exp  er iment  s  o  n  Meande ring  Channels  with  Flood  Pla ins 

Since  the  primary  objective  was  to  explore  the 
characteristics  of  the  flow  in  the  meandering  channel  when 
the  flow  tops  its  banks,  the  flume  and  the  equipment  were 
set  up  for  experiments  in  the  curved  channel  in  a  straight 
valley.  The  discharge  was  increased  until  the  depth  of  flow 
on  the  flood  plain  reached  about  one  inch  (d=2.54  cm). 
Uniform  flow  in  the  straight  reach  of  the  flume  was  ensured 
by  manipulating  the  tailgate  and  observing  the  depth  of  flow 
in  the  different  reaches. 

The  experimental  arrangement  and  the  instrumentation 
were  the  same  as  for  the  curved  channel.  While  the  velocity 
profiles  in  the  main  channel  were  measured  in  the  same  way 
as  previously  explained  in  curved  channels  without  flood 
plain,  the  profiles  on  the  flood  plain  were  spaced  0.1  feet 
(3.0  cm)  apart . 
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Two  experiments  in  the  meandering  channel  with  a  flood 
plain  were  conducted.  The  first  experiment  had  the  depth  in 
the  main  channel  equal  to  0.20  feet  (D=6.1  cm).  Depth  of  the 
flood  plain  flow  was  about  0.08  feet  (d=2.4  cm).  For  this 
run,  cross-section  numbers  1  to  15  were  investigated  in 
detail.  For  the  second  run,  where  D=0.235  feet  (7.2  cm)  and 
d=0.11  feet  (3.35  cm),  the  velocity  measurements  were 
carried  out  for  all  the  cross-sections  (numbers  1  to  15)  but 
only  at  a  level  of  y'=0.05  feet  (1.5  cm) 


5 . 5  F low  Visualization  Results 

Red  dye  with  the  same  density  as  water  was  injected 
into  the  flow  through  a  tube  of  1.5  mm  OD  and  1.0  mm  ID, 
with  a  constant  head  of  0.4  inch  (1.0  cm)  in  order  to 
clearly  observe  the  flow  path.  To  identify  what  is  happening 
at  different  levels  above  the  bed,  three  levels  were  chosen: 
one  located  directly  on  the  bed,  one  at  half  the  depth  and 
another  at  the  water  surface.  Cross-section  number  1  (at  the 
entrance  of  the  second  loop)  was  used  as  an  injection 
section.  Starting  from  the  flood  plain  bank,  the  red  dye 
spread  as  in  any  normal  test  in  a  simple  rectangular 
channel.  This  part  of  the  test  where  the  dye  spread  up  to  a 
point  of  about  9  inches  (c^=22.9  cm)  from  the  bank,  showed 
that  there  is  no  measurable  secondary  current  in  this 
region.  Consequently,  this  part  of  the  flood  plain  channel 
is  the  minimum  width  of  the  flood  plain  (c^=9  inches,  22.9 
cm),  where  the  flow  in  this  area  follows  the  alignment  of 
the  valley,  which  is  straight  (see  photos  in  Figure 
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5.10(a)). 

In  the  main  channel,  when  the  dye  was  injected  below 
the  flood  plain  level,  helicoidal  motion  was  observed  (see 
photo  in  Figure  5.10(b))  as  in  the  case  of  a  curved  channel 
without  a  flood  plain.  But,  the  helicoidal  motion  was  not 
confined  to  the  main  channel.  When  a  particle  comes  up  above 
the  flood  plain,  the  spiral  motion  will  continue  on  to  the 
flood  plain  area.  If  the  dye  is  injected  above  the  flood 
plain,  then  the  water  has  a  tendency  to  move  more  or  less 
straight  in  a  width  equal  to  the  amplitude  of  the  meander 
belt  (see  photo  in  Figure  5.10(c)).  Therefore,  it  can  be 
seen  that  the  width  of  the  main  channel  for  the  upper 
portion  is  equal  to  C2=20  inches  (50.8  cm).  The  other 
interesting  observation  was  that  the  particle  which  is 
located  just  below  the  flood  level,  instead  of  going  around 
the  bend,  simply  went  over  the  flood  plain.  Figures  5.  10  (a 
to  c)  illustrate  these  mechanisms  of  the  flow. 


5  . 6  Data  Dandling  and  Presentation  o  f  Pesul ts 

In  Chapter  II  the  data  acquisition  system  was 
described.  It  was  used  to  record  the  information  for  each 
test.  Then  the  data  processing  and  analysis  were  carried  out 
using  the  University  of  Alberta  Computer.  Due  to  some 
difficulties  in  Surface  II  Program,  in  regions  of  very  steep 
gradients,  sometimes  the  contours  cross  each  other.  Hence, 
the  contours  should  be  interpreted  with  some  care. 
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a;  Flow:  Flood  Plain  Channel 


b;  Lower  Part  of  Flow:  Helicoidal  Motion 
in  the  Meandering  Channel 


c;  Upper  Part  of  Flow:  Straight  Flow 
in  the  Meandering  Channel 

8“f*  f5:10»PI’78raph  1  Efective  z°ne  of  Momentum  Exchange 
the  Meandering  Channel  with  Flood  Plain  Channels 
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5.6.1  Velocity  Pis  t r ibu  tion 

The  angle  of  the  velocity  vector  and  the  distribution 
of  the  radial,  vector  and  tangential  velocities  within  the 
selected  channel  cross-sections  (numbers  1  to  15)  for  run 
number  1  are  presented  in  Figures  5.11  to  5.14.  Since  every 
cross-section  had  at  least  40  vertical  velocity  profiles,  it 
was  decided  to  select  a  few  representive  profiles  for  the 
sake  of  brevity. 

A  quick  glance  at  the  vertical  distribution  of  9?  the 
angle  of  the  velocity  with  the  tangential  direction  on  the 
flood  plain  for  the  cross-section  1  to  15  in  Figures  5.11(1) 
to  5.11(15)  reveals  that  9  ranges  between  -2  to  +2  degrees 
except  for  the  region  close  to  the  meandering  main  channel. 
The  precision  of  the  yaw  probe  in  measuring  the  angle  of  the 
flow  direction  for  9  close  to  zero  degree  could  be  a  few 
degrees  (Rajaratnam  &  Muralidhar,  1969).  Using  the  modified 
procedure  of  the  yaw  probe  technique,  it  is  estimated  this 
inaccuracy  could  be  reduced  to  a  couple  degrees.  Considering 
the  above  inaccuracy,  it  is  concluded  that  the  velocity 
vector  in  the  flood  plain  channel  is  parallel  to  the  x  axis. 

Since  the  lower  part  of  the  flow  in  the  meandering 
channel  is  bounded  by  the  side  walls,  the  flow  is  forced  to 
flow  more  or  less  along  the  stream -lines  of  the  curved 
channel.  In  the  upper  part  of  the  flow,  because  of  the 
absence  of  side  walls,  the  flow  has  a  tendency  to  move 
parallel  to  the  x  direction.  The  yaw  probe  in  the  meandering 
channel  was  always  kept  parallel  to  the  local  tangential 
direction  of  the  curvature  of  the  main  ch  ann  el  at  each 
section.  Consider  the  profiles  located  at  one  inch  (2.54  cm) 
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Figure  5.1 1(1)  Angle  of  Velocity  Vector  Profiles:  Cross-Section  No. 
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Figure  5.1 1(5)  Angle  of  Velocity  Vector  Profiles:  Cross-Section  No. 
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Figure  5.11(6)  Angle  of  Velocity  Vector  Profiles:  Cross-Section  No. 
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Figure  5.1 1(7)  Angle  of  Velocity  Vector  Profiles:  Cross-Section  No. 
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Figure  5.1 1  (8)  Angle  of  Velocity  Vector  Profiles:  Cross-Section  No. 
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Figure  5.11(9)  Angle  of  Velocity  Vector 
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Figure  5.11(10)  Angle  of  Velocity  Vector  Profiles:  Cross-Section  No.  10 
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Figure  5.11(11)  Angle  of  Velocity  Vector  Profiles:  Cross-Section  No.  11 
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Figure  5.1 1  (13)  Angle  of  Velocity  Vector  Profiles:  Cross-Section  No.  13 
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Figure  5.11(14)  Angle  of  Velocity  Vector  Profiles:  Cross-Section  No.  14 
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Figure  5.11(15)  Angle  of  Velocity  Vector  Profiles:  Cross-Section  No.  15 
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Figure  5.12(1)  Radial  Velocity  Profiles:  Cross-Section  No. 
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Figure  5.12(3)  Radial  Velocity  Profiles:  Cross-Section  No. 
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Figure  5.12(5)  Radial  Velocity  Profiles:  Cross-Section  No. 
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Figure  5.12(6)  Radial  Velocity  Profiles:  Cross-Section 
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Figure  5.12(7)  Radial  Velocity  Profiles:  Cross-Section  No. 
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Figure  5.12(8)  Radial  Velocity  Profiles:  Cross-Section  No. 
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Figure  5.12(9)  Radial  Velocity  Profiles:  Cross-Section  No. 
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Figure  5.12(10)  Radial  Velocity  Profiles:  Cross-Section  No.  10 
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Figure  5.12(13)  Radial  Velocity  Profiles:  Cross-Section  No.  13 
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Figure  5.12(14)  Radial  Velocity  Profiles:  Cross-Section  No.  14 
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Figure  5.12(15)  Radial  Velocity  Profiles:  Cross-Section  No.  15 
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Figure  5.13(1)  Velocity  Vector  Profiles:  Cross-Section  No. 
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Figure  5.13(3)  Velocity  Vector  Profiles:  Cross-Section  No. 
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Figure  5.13(6)  Velocity  Vector  Profiles:  Cross-Section  No. 
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Figure  5.13(7)  Velocity  Vector  Profiles:  Cross-Section  No. 
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Figure  5.13(8)  Velocity  Vector  Profiles:  Cross-Section  No. 
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Figure  5.13(9)  Velocity  Vector  Profiles:  Cross-Section  No. 
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Figure  5.13(10)  Velocity  Vector  Profiles:  Cross-Section  No.  10 
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Figure  5.13(1 1)  Velocity  Vector  Profiles:  Cross-Section  No.  11 
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Figure  5.13(13)  Velocity  Vector  Profiles:  Cross-Section  No.  13 
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Figure  5.13(14)  Velocity  Vector  Profiles:  Cross-Section  No.  14 
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Figure  5.13(15)  Velocity  Vector  Profiles:  Cross-Section  No.  15 
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Figure  5.14(3)  Tangential  Velocity  Profiles:  Cross-Section  No. 
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Figure  5.14(6)  Tangential  Velocity  Profiles:  Cross-Section  No. 
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Figure  5.14(7)  Tangential  Velocity  Profiles:  Cross-Section  No. 
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Figure  5.14(10)  Tangential  Velocity  Profiles:  Cross-Section  No.  10 
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Figure  5.14(1 1 )  Tangential  Velocity  Profiles:  Cross-Section  No.  11 
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Figure  5.14(13)  Tangential  Velocity  Profiles:  Cross-Section  No.  13 
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Figure  5.14(14)  Tangential  Velocity  Profiles:  Cross-Section  No.  14 
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Figure  5.14(15)  Tangential  Velocity  Profiles:  Cross-Section  No.  15 


169 


from  the  side  wall  in  the  main  channel  for  cross-sections 
number  1  to  15  .  One  could  clearly  see  from  the  profiles  in 

Figures  5.11(1)  to  5.11(15)  that  the  angle  of  the  velocity 
vector  with  respect  to  the  local  tengential  direction  9, 
varies  a  little  for  the  lower  part  of  the  flow  (y=0.0  to  h), 
but  for  the  upper  part  of  the  flow  (y=h  to  water  surface), 

0  varies  in  such  a  way  that  the  velocity  vector  is  parallel 
to  x  axis.  For  example,  for  FI  in  Figure  5.11(3),  where  the 
cross-section  is  located  at  30  degrees  with  respect  to  z,  9 
has  value  of  a  few  degrees  when  y=0.0  up  to  0.1  feet  and 
suddenly  9  increases  to  a  constant  value  of  about  -30 

degrees  from  y=h=0.125  feet  to  P=0.20  feet. 

It  can  be  seen  from  Figure  5.12  that  the  component  of 
the  velocity  vectors  in  z  direction  (w)  in  the  flood  plain 
channel  are  negligible  except  in  the  region  near  the  main 
channel  which  is  affected  by  the  secondary  current  in  the 

meandering  channel. 

Since  the  upper  part  of  the  flow  in  the  meandering 
channel  is  flowing  almost  parallel  to  the  x  axis,  the 
magnitude  of  vr  is  very  small  whereas  in  the  lower  part  of 
the  main  flow,  the  magnitude  of  vr  reaches  up  to  40  percent 

of  the  velocity  vector. 

The  tangential  and  vector  velocity  profiles  in  the 
flood  plain  channel  were  replotted  with  u/u*  versus 
l0iy  "*-  ,  where  u  simply  replaced  V  and  u0  because  there 
were  no  difference  between  the  two  velocities.  From  these 
plots  (these  are  not  presented  here),  it  was  found  that  the 
velocity  distribution  agreed  well  with  the  well  known 
Karman-Prandtl  eq ua t ion .  ( Sch 1 i ch t ing ,  1  968)  for  a  smooth 
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even  for  the  region  near  the  main  channel.  It  was  noticed 
that  near  the  free-surface  there  was  a  dip  in  the  velocity 

profile  (Fajaratnam  &  Muralidhar,  1968). 

Figures  5.13  and  5.14  show  that  V  and  increase  with 

y,  the  distance  from  the  bed  of  the  meandering  channel 
rather  rapidly  and  then  assume  an  almost  constant  value  up 
to  the  flood  plain  level.  From  this  point  on,  if  the  profile 
is  close  to  the  outermost  wall,  the  velocity  decreases 
continuously  to  eventually  become  approximately  constant 
near  the  water  surface.  Likewise,  if  the  profile  is  near  the 
innermost  wall,  the  velocity  increases  continuously  to 
become  constant  near  the  water  surface.  This  behavior  was 
noticed  in  the  straight  main  channel  with  straight  flood 

plains  (Chapter  III  and  IV). 

From  the  above  figures  it  was  concluded  that  separating 

the  effect  of  the  momentum,  exchange  due  to  the  compound 
cross-section  fro,  the  one  due  to  the  effect  of  the  curved 
channel  is  difficult,  if  possible  at  all.  A  comparison  of 
U(y)  for  flow  with  and  without  the  flood  plain  shows  a 
little  difference  for  the  lower  part  of  the  profile.  For  the 
upper  part  ,  above  the  flood  plain  level,  as  we  should 
expect  from  the  experience  in  straight  channel,  it  should 
have  a  dip  with  approximately  a  constant  velocity  for  each 
vertical  profile.  The  dips  in  meander  flows  are  not  similar 
to  the  dips  in  the  straight  compound  channels. 


Figure  5.15  shows  the  lateral  distribution  of  the 
tangential  velocity  (u(z))  for  run  number  1  at  different 
levels  above  the  flood  plain  (y ' ) -  For  run  number  2,  u&(z) 
is  presented  for  y'=0.05  feet  (1.5  cm)  is  shown  in  Figure 
5.16.  From  Figure  5.15  it  can  be  seen  that  ug  increases  with 
z,  the  distance  from  the  side  wall  of  the  flood  plain  until 
it  reaches  almost  a  constant  value.  Eue  to  the  interaction 
between  the  curved  channel  and  the  flood  plain,  u@  increases 
to  a  maximum  value  at  the  innermost  junction  plane  between 
the  main  channel  and  flood  plain.  The  tangential  velocity  in 
the  main  channel  decreases  rather  rapidly  toward  the 
outermost  junction  plane  between  the  meandering  channel  and 
flood  plain  channel.  At  the  outermost  junction  plane,  one 
could  see  a  small  increase  in  uQ  in  the  interaction  zone, 
but  eventually  u@  approaches  a  constant  value  and  then  will 
again  decrease  as  the  bank  of  the  flood  plain  is  approached. 
Due  to  the  uniformity  of  the  velocity  profile  in  the  upper 


part  of  the  flow  (u&(y)>,  the  curves  of  the  lateral 
distribution  of  u@  in  the  main  channel  for  different  levels 

above  the  flood  plain  level  are  the  same. 

Figures  5.15  and  5.16  clearly  show  that  the  velocity 


near  the  edge  of  the  main  channel  and  the  flood  plain  has 
been  increased  because  of  the  effect  of  the  interaction 
between  the  two  channels.  As  can  be  seen,  the  lateral 
distribution  of  the  velocity  over  the  flood  plain  on  both 
sides  of  the  main  channel  is  not  the  same.  On  the  convex 
side,  over  the  flood  plain,  the  velocity  steadily  increases 
up  tl  the  crest  of  the  meander.  Then  it  decreases  until  it 
reaches  the  meandering  channel.  However,  the  flow  velocity 
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Figure  5.15  Continued 
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Figure  5.16  Lateral  Tangential  Velocity  Distribution:  Cross-Section  No.  1  to  15  (RUN2) 
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over  the  flood  plain  on  the  convex  side  is  always  higher 
than  on  the  concave  side.  Therefore,  it  was  decided  to  draw 

isovel  waps  for  different  levels. 

From  Figure  5.11,  which  reveals  the  distribution  of  9, 

one  can  conclude  from  the  sudden  change  of  9  at  y=h,  that 
the  velocity  vector  direction  is  changing  from  its  direction 
in  the  main  channel  to  becoming  almost  parallel  to  the  x 
axis  at  the  flood  plain  elevation.  The  sudden  change  of  the 
direction  and  the  magnitude  of  V  suggests  there  is  a 
rotational  shear  as  well  as  a  parallel  shear  in  the 
horizontal  (parallel  to  the  meandering  bed)  plane. 

5.6.2  Cross-Sect ional  Isovel 

For  cross-section  numbers  1  to  15,  maps  of  isovels  for 
the  tangential  velocity  were  prepared.  These  cross-sectional 
isovel  maps  are  shown  in  Figure  5.17.  These  maps  show  two 

important  mechanisms  of  the  flow. 

Firstly,  in  the  main  channel  below  the  flood  plain 

level,  the  higher  velocity  will  occur  near  the  convex  side 
(the  innermost  wall)  at  the  beginning  of  the  bend  and  will 
shift  toward  the  concave  (outermost  wall)  with  increasing  9. 
Nevertheless,  it  will  never  occur  close  to  the  concave  wall. 
The  amount  of  shifting  of  u0  can  be  seen  by  comparing  the 
cross-section  number  5  with  cross-section  number  11.  They 
are  located  45  degrees  before  the  crest  and  past  the  crest 

of  the  meandering  channel  respectively. 

Secondly,  the  meandering  channel  flow  will  affect  the 
flood  plain  channel  flow.  As  the  result  of  the  lateral 
momentum  exchange  between  the  two  channels,  increased 
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Figure  5.17  Continued 
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Figure  5.17  Continue 
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velocities  will  occur  on  the  flood  plain  area.  The 
approximate  zone  of  interaction  can  be  distinguished  in  the 
channel.  Once  again,  comparison  of  the  two  cross-sections 
(numbers  5  and  11)  will  show  that  they  are  not  identical. 

They  differ  mainly  with  regards  to  the  acceleration  and 
deceleration  of  the  flow  on  the  flood  plain  and  the  effect 
of  9.  The  isovel  maps  for  the  main  channel  at  these  sections 

are  not  the  same. 

5.6.3  I s  ove 1  Nap 

For  the  entire  second  loop  at  different  levels  above 
the  flood  plain,  u0,  *,  V  and  a),  (where  *  the  angle  of  the 

velocity  vector  with  the  x  axis),  contours  are  shown  in 
Figures  5.18  ,  5.19,  5.20  and  5.22  for  run  number  1.  For  run 
number  2  the  same  contour  maps  for  the  level  y'=0.05  feet 

(1.15  cm)  are  shown  in  Figures  ^.21,  5.22. 

Figures  5.19,  and  5.21  show  that  the  deviation  of  the 
angle  of  the  vector  velocity  from  the  x  axis  ranges  from.  -20 
to  +20  degrees,  which  mostly  covers  the  main  channel  and  its 
surrounding  area.  These  maps,  by  themselves,  suggest  that 
the  radial  velocity  component  in  the  main  channel  and  its 
neighbourhood  in  the  flood  plain  (interaction  zone)  are 
large  by  comparison  with  the  flow  on  the  flood  plain  area 

far  away  from  the  main  channel. 

Figures  5.19  and  5.21  also  show  the  variation  of  the 

angle  of  deviation  of  the  velocity  vector  *,  from  the  x 
direction,  (*  being  negative  if  the  rotation  of  V  to  x  is 
clockwise).  Figure  5.22  shows  isovel  maps  of  the  radial 
velocity  component.  These  contour  maps  are  given  for  y' 
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=0.003  and  0.04  feet  (0.9  and  12.2  mm)  for  run  number  1  and 
y ' =  0 . 0  5  feet  (1.52  cm)  for  run  number  2.  Figures  5.19  and 
5.22  provide  evidence  of  the  flow  exchanges  between  the 
meandering  main  channel  and  the  flood  plain,  confirming  the 
existence  of  secondary  currents  and  showing  that  the 
magnitude  of  the  lateral  velocity  component  could  reach  (at 
some  locations)  be  as  much  as  50%  of  V.  One  should  also 
notice  the  divergent  and  convergent  patterns,  even  m  the 
straight  channel. 

Tangential  and  vector  velocity  contour  maps  looked 
alike  for  each  y'  and  therefore,  V  isovels  are  presented  in 
Figure  5.20  for  only  three  levels  y'=0.002,  0.04  and  0.07 

feet  (y ' =  0  •  6  mm,  1.2  cm  and  2.2  cm)  The  isovel  maps  in 
Figure  5.20  and  5.22  show  that  the  maximum  velocity  contour 
is  located  near  the  convex  wall,  but  they  are  not 
symmetrical  with  respect  to  the  crest  of  the  meandering 
curve.  This  is  because  the  behaviour  of  the  flow  approaching 
the  crest  is  not  the  same  as  the  behaviour  of  the  flow  past 
the  crest.  Another  comment  can  be  made  about  the 
acceleration  and  deceleration  of  the  velocity  vector  by 
noticing  the  closed  contours  for  equal  values  of  (J»  . 

Figure  5.22  shows  contours  of  the  lateral  velocity 
component  at  level  y'=0.003  and  0.04  feet  (0.9  mm  and  12.2 
mm)  for  ran  number  1  and  y'=0.05  feet  (15.2  mm)  for  run 
number  2.  This  figure  clearly  shows  the  zones  of  the  flood 
plain  disturbed  by  the  main  meandering  channel  and  also  the 

magnitude  of  the  secondary  current. 

It  was  concluded  that  the  helicoidal  currents  were  not 

channel  and  spread  to  a  broader  area 


confined  to  the  main 


. 
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equal  to  the  width  of  the  meander  loop.  The  exact 
contribution  of  the  momentum  exchange  from  the  curved 
channel  and  the  compound  cross-section  to  the  lateral 
velocity  redistribution  is  not  known.  But,  qualitatively, 
the  effect  of  the  curved  channel  and  the  exchange  of 
momentum  for  u(z)  could  be  seen  from  lateral  tangential 

velocity  profiles. 


5  .  7  Discussion  on  Interaction  Results 

There  is  an  absence  of  sound  knowledge  about  the  flow 
in  curved  channels,  even  in  the  simplest  cross-section 
(rectangular).  Addition  of  flood  plains  to  meandering 
channels  brings  additional  difficulties.  Tn  light  of  the 
study  described  in  chapters  III  and  IV,  transfer  of  momentum 
will  take  place  from  the  region  of  high  velocity  to  the 
region  of  the  lower  velocity  (in  the  case  of  straight 
channel  from  the  main  channel  to  the  flood  plain  channel). 

In  the  case  of  meandering  channels  with  flood  plains,  one 
should  expect  the  transfer  of  momentum  to  occur  between  the 
high  and  low  region  of  the  velocity. 


The  secondary  current,  which  is  generated  in  the 
meandering  channel  because  of  the  curvature  of  the  main 
channel,  will  transport  momentum  in  and  out  from  the 
meandering  channel  to  the  flood  plain.  The  helocoidal  motion 
in  the  lower  part  of  the  curved  channel  is  confined  in  the 
main  channel  because  of  the  presence  of  the  solid  walls. 
However,  in  the  upper  part  of  the  flow,  the  absence  of  solid 
boundaries  allows  the  helicoidal  motion  to  continue  on  to 
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the  flood  plain.  Therefore,  in  Figure  5.22  an  isovel  man  for 
the  lateral  component  of  the  velocity  vector  vr>  f°r  a  few 
selected  cross-sections  is  shown.  The  secondary  currents  are 
mostly  confined  to  a  region  equal  to  the  total  width  of  the 
meander . 

The  flow  above  the  flood  plain  level  has  a  tendency  to 
flow  more  or  less  in  the  x  direction.  Because  of  changes  m 
the  depth  of  the  flow  in  x  direction,  the  flow  will 
accelerate  in  the  deeper  part  Con  the  main  channel)  and  will 
slowly  decelerate  on  the  shallow  part  (in  the  flood  plain) 
by  excess  friction.  Comparison  of  the  lateral  tangential 
velocities  of  cross-section  numbers  1  to  15  with  the 
symmetrical  cross-section  number  0  (i.e.  see  Figure  5.16) 

will  bring  the  following  points  to  light: 

1.  increased  velocity  in  the  flood  plain  channel  due  to 
exchange  of  momentum  between  the  main  flow  and  the  flood 

plain  f low , 

2.  decreased  velocity  in  the  main  channel  due  to 
interaction.  However,  the  lateral  tangential  velocities 
retain  their  features  such  as  higher  velocities  near  the 
innermost  wall  and  lower  near  the  outermost  wall, 

3.  due  to  the  acceleration  and  the  deceleration  mechanisms, 
at  the  beginning  of  the  bend,  the  velocity  on  the  flood 
plain  will  increase  up  to  cross-section  number  3  (the 
crest  cross-section)  and  decrease  gradually, 

4.  the  sudden  change  of  the  direction  of  the  velocity 
velocity  in  the  main  channel  at  plane  z=h  introduces  a 

rotational  shear  in  this  plane. 

In  most  cross-sections  the  lateral  velocity  profiles  near 
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the  convex  edge  of  the  flood  plain  bears  the  maximum 
velocity  which  normally  is  associated  with  high  shear 
stress.  Therefore,  at  high  flood  flow  the  convex  edge  of  the 
flood  plain  will  be  eroded  and  the  channel  will  straighten 

out. 


5  .  p,  Conclusions  and  Comment  s  on  Further  Direct  ion. 

The  present  study  of  flow  field  in  a  smooth  rigid 
boundary  model  of  a  curved  channel  and  a  meandering  channel 
in  a  straight  flood  plain  valley  has  led  to  the  following 
qualitative  conclusions. 

The  logarithmic  velocity  distribution  law  may  be 
applicable  only  in  the  immediate  region  close  to  the  bed  of 
the  curved  channels  (about  20 7,  of  the  depth  of  flow).  It  is 
shown  that  the  distribution  of  uQ(y)  in  the  core  of  flow  is 
mostly  uniform.  The  magnitude  of  u@  is  dependent  on  Cos6. 
Also,  there  is  a  relationship  between  P,  u@  and  6.  The 
normalized  lateral  velocity  distribution  u@(z),  for 
different  cross-sections  has  been  studied.  They  are  similar 
if  the  distribution  curves  rotate  around  the  point  of 
intersection  of  all  distribution  curves,  where  it  is  located 
at  20%  of  B  from  the  innermost  wall  of  the  curved  channel. 

The  shear  stress  near  the  convex  region  of  the 
meandering  channel  was  about  2.5  times  as  large  as  that  for 
a  uniform  channel  of  the  same  width  and  for  the  same 

hydraulic  radius  and  discharge. 

The  interaction  between  the  inbank  and  overbank  flows 

results  in  lateral  momentum  transfer  which,  in  turn. 
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produces  increased  flood  plain  velocity  and  increased  bed 
shear  stresses.  This  interaction  also  causes  reductions  in 
the  meandering  channel  velocity.  The  lower  part  of  the  flow 
in  the  main  channel  bounded  by  the  side  walls,  flows  more  or 
less  parallel  to  the  boundary.  The  upper  part  of  the  flow 
above  the  flood  plain  level  flows  parallel  to  the  x 
direction.  This  change  of  the  direction  of  the  flow  in  a 
horizontal  plain  will  cause  an  additional  rotational  shear 
between  the  adjacent  horizontal  plains.  At  this  stage  of  the 
study,  no  attempt  was  made  to  determine  the  magnitude  of  the 
rotational  shear  stresses  nor  the  effect  on  reduction  of  the 

discharge  due  to  this  extra  shear. 

An  interesting  observation  was  the  acceleration  of  the 

flow  from  the  main  channel  to  the  flood  plain  up  to  the 
middle  of  the  loop  in  the  flood  plain  area.  At  this  point 
the  flow  started  to  decelerate  until  it  reached  the  end  of 

the  loop. 


5  .  9  Further  Direction 

The  existing  models  for  prediction  of  the  flow  in  a 
curved  open  channel  are  not  adequate  for  describing  the 
flow.  This  is  the  case  mainly  because  the  models  ignore  the 
fact  that  9  is  one  of  the  most  important  parameters.  Besides 
that,  the  assumptions  made  to  simplify  the  equation  of 
motion  can  not  be  justified  (fully  developed  flow, 
logarithmic  law,  etc.). 

Vith  this  state  of  knowledge  for  flow  in  curved 
channels,  the  addition  of  a  flood  plain  makes  the  problem 
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even  more  difficult.  Pence  it  might  be  advisable  to  first 
tackle  the  problem  of  flow  in  curved  channels  and  then 
attempt  to  solve  the  problem  of  flow  in  curved  main  channels 

with  flood  plains. 


°'3 
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6.  Summary  of  Conclusions 


6  .  l  Concep  t  o  f  Interact  ion  between  Main  Ch  anneJL  and  Flo°  j. 
Plain  Flows 

The  concept  of  the  momentum  transfer  from  high  to  low 
velocity  region  by  turbulent  mixing  was  tested  in  the  case 
of  a  main  channel  (region  of  high  momentum)  with  a  straight 
flood  plain  (region  of  low  momentum).  The  lateral  transfer 
of  momentum  from  the  main  channel  to  the  flood  plain  channel 
results  in  the  creation  of  a  mixing  region  between  the  mam 

channel  and  the  flood  plain. 

The  velocity  and  boundary  shear  stress  distributions 

were  measured  for  a  range  of  flow  conditions  in  a 
rectangular  compound  cross-section  channel  with  both 
straight  symmetrical  and  asymmetrical  flood  plains  and  also 
for  a  meandering  main  channel  with  straight  flood  plains.  A 
series  of  figures  show  the  relationships  between  boundary 
shear  stress,  aspect  ratio  and  relative  roughness  for 
rectangular  channels.  Methods  of  estimating  the  local 
boundary  shear  stress,  velocity  and  discharge  for  a  simple 
rectangular  channel  and  a  compound  cross-section  channel 

shape  are  described. 


6^2  Vide  Main  Channel  and  Vide  Flood  Plain  Channel 

For  a  straight  wide  main  channel  with  a  wide  flood 
plain,  the  interaction  between  the  main  channel  and  the 
flood  plain  flows  has  been  studied  experimentally.  The 
results  have  demonstrated  the  transport  of  the  longitudinal 
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momentum  from  the  main  channel  to  the  flood  plain.  The 
lateral  velocity  profiles  in  the  flood  plain  and  the  mam 
channel  above  the  level  of  the  flood  plain,  if  viewed  with 
respect  to  the  undisturbed  flood  plain  velocity,  have  been 
found  to  be  approximately  similar.  Empirical  correlations 
have  been  found  for  the  velocity  and  length  scales.  It  seems 
the  structure  of  the  flow  in  a  main  channel  with  a  flood 
plain  could  be  analysed  using  the  concept  of  lateral 
momentum  transfer  and  by  treating  the  interaction  region  as 
a  turbulent  shear  layer.  The  vertical  velocity  profiles  can 
be  described  by  logarithmic  velocity  distribution  law.  This 
suggests  analytical  possibilities  of  predicting  the 
magnitude  of  the  turbulent  eddy  viscosity  for  the 
interaction  region. 


The  bed  shear  in  the  flood  plain  near  the  junction 
plane  with  the  main  channel  is  increased  considerably  and 
the  main  channel  bed  shear  is  decreased  because  of  the 
interaction  between  the  main  channel  and  the  flood  plain 
flows.  These  changes  have  been  analysed  using  similarity 
ideas.  Some  estimates  of  this  penetration  distance  could 
perhaps  be  made  using  the  results  on  turbulent  compound 
shear  layers.  The  scale  factors  for  the  shear  stress 
distribution  have  been  developed  and  D/d  has  been  found  to 


,e  the  main  parameter  to  indicate  the  extent  of  the 
interaction  between  the  main  channel  and  flood  plain  flows. 

There  is  a  reduction  in  carrying  capacity  in  the 
compound  cross-section  channel  due  to  the  flow  interaction. 
It  is  possible  to  express  this  reduction  of  the  discharge 
analytically.  The  results  of  one  experiment  with  a  rough 
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flood  plain  is  included. 


6  .  3  Meandering  Channel  in  a  S  t  ra  igh  t  Va  Hey 

Flow  in  a  curved  channel  and  meandering  channel  in  a 
straight  flood  plain  valley  was  also  studied. 

The  logarithmic  velocity  distribution  law  may  be 
applicable  only  in  the  immediate  region  close  to  the  bed  of 
the  curved  channel  for  about  20%  of  the  depth  of  flow.  It  is 
shown  that  the  vertical  distribution  of  the  velocity  in  the 
core  of  flow  is  mostly  uniform.  The  magnitude  of  the 
velocity  is  dependent  on  the  angular  position  of  the 
section.  Normalized  lateral  velocity  distributions  for 
different  cross-sections  are  similar.  The  bed  shear  stress 
near  the  convex  region  of  the  meandering  channel  was  about 
to  2.5  times  as  large  as  that  for  a  uniform  channel  of  the 
same  width  and  for  the  same  hydraulic  radius,  roughness  and 

discharge  . 

The  interaction  between  the  inbank  and  overbank  flows 
results  in  lateral  momentum  transfer  which,  in  turn, 
produces  increased  flood  plain  velocity  and  increased  the 
bed  shear  stresses.  This  interaction  also  causes  reductions 
in  the  meandering  channel  velocity  and  bed  shear  stresses. 
The  lower  part  of  the  flow  in  the  main  channel  bounded  by 
the  side  walls  flows  more  or  less  parallel  to  the  boundary. 
The  upper  part  of  the  flow  above  the  flood  plain  level  flows 
parallel  to  the  x  direction.  This  change  of  the  direction  of 
the  flow  in  horizontal  plain  will  cause  an  additional 
rotational  shear  between  the  adjacent  horizontal  planes.  At 
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this  stage  of  the  study,  no  attempt  was  made  to  determine 
either  the  magnitude  of  the  rotational  shear  stresses  or  the 
effect  on  reduction  of  the  discharge  due  to  this  extra 
shear  . 


The  flow  accelerates  from  the  main  channel  to  the  flood 
plain  up  to  the  middle  of  the  loop  in  the  flood  plain  area. 
At  this  point  the  flow  starts  to  decelerate  until  it  reaches 
the  end  of  the  loop.  The  present  results  on  the  meandering 
channel  with  flood  plains  could  be  considered  to  constitute 
only  an  exploratory  study. 
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NOTATIONS 

The  following  symbols  are  used  in  this  report: 


B  , 


C, 


D 

d 

d' 


F 

f 

g 

h 


h' 


h 

h 


o 

1 


k  and  k,  to  k,  =  constants 

1  6 

=  width  of  the  main  channel 
=  width  of  flood  plain  channel 
=  length  scale  for  flood  plain 

=  length  scale  for  main  channel 

=  total  length  scale 

=  shear  length  scale 

=  width  of  first  flood  plain  channel 
=  width  of  second  flood  plain  channel 
=  minimum  width  of  the  flood  plain  channel 

=  maximum  width  of  the  flood  plain  channel 

=  depth  of  main  flow 
=  depth  of  flow  in  flood  plain 

=  depth  of  flow  in  flood  plain  at  second  level 

in  multi-stage  flood  plain 
=  Froude  number 
=  function 

=  acceleration  due  to  gravity 

=  depth  of  main  channel,  below  the  level  of 
flood  plain 

=  depth  of  main  channel  below  the  second  level 
of  the  flood  plain  in  mult i-s tage f lood  plain 

=  static  pressure  head 

=  dynamic  pressure  head  of  hole  No.l  of  the 
y  aw  probe 
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dynamic  pressure  head  of  hole  No. 2  of  the 
yaw  probe 

dynamic  pressure  head  of  hole  No. 3  of  the 
yaw  probe 

Nikuradse  sand  roughness 

value  of  z'  where  At  =0.0 

o 

length  of  meander  channel 

entrance  length 

exit  length 

wave  length 

pressure 

static  pressure 

pressure  of  hole  No.l  of  the  yaw 
probe 

pressure  of  hole  No .2  of  the  yaw 
probe 

pressure  of  hole  No. 3  of  the  yaw 

probe 

discharge 

centerline  radius  of  bend 
innermost  radius  of  bend 
outermost  radius  of  bend 
slope 

slope  of  channel  (bed  slope) 

water  surface  slope  in  radial  direction 

water  surface  slope 

length  of  tangent 

reference  velocity 

velocity  in  the  longitudinal  (x)  direction 
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u 

m 

* 

u 

m 

u 

o 

u 

oo 


u  . 

*oo 

V 

v  o  r  w 
r 

x 


z 


z 


* 


n  ,  o  ' 

© 

9 

V 

k 


om 

T' 


=  tangential  velocity 
=  value  of  u  for  z  =  0 . 0 
=  value  of  u  for  z=B/2 
=  value  of  u  for  y=h 

=  undisturbed  value  of  u  in  the  flood  plain 
for  any  y  ' 

=  shear  velocity 

=  undisturbed  shear  velocity 

=  velocity  vector 

=  radial  velocity 

=  longitudinal  distance 

=  vertical  distance  from  bed  of  main  channel 

=  vertical  distance  from  bed  of  flood  plain 

=  loction  of  the  dip  in  the  main  channel 
=  lateral  distance  from  the  centerline  of  main 
channel 

=  lateral  distance  for  flood  plain 
(  z  =  B  /  2  ,  z  '  =  0 . 0  ) 

=  dimensionless  lateral  distances 
=  angular  coordinate 

=  angle  of  attack  for  the  yaw  probe 

=  central  angle  of  bend 

=  kinematic  viscosity  of  water 

=  Von  Kantian  universal  constant 

=  bed  shear  stress 

=  T  at  z  =  C . 0 
o 


ora 


at  z  =  B  /  2  or  z ' =0 . 0 


^  \ 

. 
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T  0  GO 

=  undisturbed  TQ  in  the  flood  plain 

> 

H 

O 

=  t Q  in  flood  plain  minus 

Tom 

=  t  with  side  walls  of  the  main  channel 

om 

extended  to  total  depth 

T* 

=  average  longitudinal  shear  stress  in 

the  junction  plane 

4» 

=  angle  of  deviation  of  velocity  vector  from 

longitudinal  direction  of  channel  projected 
on  horizontal  plane,  positive  outward 

( - ) ave 

=  average  value  for  section 

D 

=  mean  value 

A 

=  wave  length 

£ 

=  eddy  viscosity 

P 

=  density  of  water 

Y 

=  specific  weight  of  water 

. 
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APPENDIX  A 


A  - 1  Turbulent  Flow  in  Pect  angular  Channels 

Perhaps  the  most  simple  cross-section  in  an  open 
channel  is  a  rectangular  one.  Even  in  a  simplified  case  like 
a  straight  channel  with  rectangular  cross-section,  the 
equations  of  motion  for  turbulent  flow,  along  with  the 
continuity  equation  cannot  be  solved.  In  order  to  formulate 
an  approximate  solution,  it  is  required  to  adopt  empirical 
assumptions.  It  is  necessary  to  supplement  the  original 
hypothesis  with  additional  hypotheses  which  may  vary  from 
case  to  case.  The  form  of  certain  functions,  or  at  least 
certain  numerical  values,  must  be  extracted  from 
experiments.  In  order  to  ignore  the  effect  of  the  side 
walls,  one  must  always  assume  the  channel  is  very  wide. 

Gruff  (1965),  and  subsequently  Rajaratnam  &  Muralidhar 
(1969),  obtained  considerable  experimental  data  for  smooth 
rectangular  channels.  Rajaratnam  &  Muralidhar  (1969),  in 
their  paper,  presented  the  variation  of  maximum  bed  shear 
stress  with  respect  to  aspect  ratio.  In  addition,  they 
expanded  upon  the  distribution  of  bed  shear  stresses  and 

wall  shear  stresses. 

This  study  presents  some  supplementary  measurements  on 
the  shear  stress  and  velocity  distribution  in  smooth  and 
rough  channels.  Together  with  the  data  available  from  Cruff 
(1965).  Rajaratnam  6  Muralidhar  (1969)  and  Ghosh  6  Jena 
(1971)  and  supplementary  measurements  in  a  smooth  and  rough 
rectangular  channels,  the  author  was  able  to  relate  the 
effect  of  wall  (aspect  ratio)  and  the  bed  roughness  with 
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aximum  bed  shear  stresses.  Furthermore,  the  distribution  of 
bed  and  wall  shear  stresses  with  width  and  depth  of  flow 
respectively  was  obtained.  Further,  the  average 
cross-sectional  velocities  for  the  entire  section  in  b 
cases,  smooth  and  rough,  were  compared  with  average 
velocities  obtained  by  dividing  discharge  by  cross-sectional 

area. 


A  _2  Experimental  Arrangements  an_d  Experiments 

The  experiments  were  performed  in  an  eight  inch  (20.2 
cm)  wide,  seven  inch  (17.8  cm)  deep  and  60  feet  (18.3  m) 
long  channel  with  adjustable  slope.  Water  entered  the  flume 
from  a  constant  head  reservoir  provided  with  suitable 
screens.  The  discharge  was  measured  by  means  of  an 
orifice-meter  located  in  the  supply  line  to  the  head-tank. 

The  height  of  the  water  surface  was  measured  with  point 
gauges.  The  water  surface  slope  was  measured  by  a  pressure 
transducer.  In  order  to  obtain  uniform  flow,  the  tail  gate 
was  adjusted  until  the  water  surface  slope  for  the  middle  20 
feet  (6.1  m)  of  the  flume,  coincided  with  the  bed  slope.  The 
test  section  was  located  in  the  middle  of  the  flume. 
Velocities  were  measured  using  a  1/8  inch  (3mm)  external 
diameter  Prandtl-tyPe  pitot  static  tube.  The  same  Prandtl 
tube  also  served  as  the  Preston  tube  for  the  smooth 

b  ound  a  ry  • 

In  the  case  of  the  smooth  boundary,  most  data  were 
taken  from  the  work  of  F.ajaratnam  &  «uralidhar  (  1  969  ). 
Complementary  runs  were  conducted  to  verify  the  average 
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curve  and  to  fill  the  gap  between  data  points. 

In  the  case  of  the  rough  bed  with  smooth  walls,  data 
were  obtained  from  a  Pepartmental  report  by  Rajaratnam 
(  1  967-1  969),  from  Ghosh  Uena  (1971)  and  from  several 
experiments  conducted  by  others.  Futher  experiments  were 
conducted  with  pre-molded  rubber  mats  with  hemispheres  of  1 
inch  (2.54  cm)  and  1/2  inch  (12.7  mm)  in  diameter. 

All  data  were  chosen  in  such  a  way  that  the  Froude 
number  varied  from  0.10  to  0.70.  Velocity  and  the  boundary 
shear  stress  distributions  were  obtained  in  one  half  of  the 
channel.  If  y  is  the  normal  distance  above  the  bed  and  if  z 
denotes  the  transverse  distance  from  the  centerplane  of  the 
channel,  then  X-z/(B/2)  and  n=y/D,  where  B  is  the  width  of 
the  channel  and  D  is  the  depth  of  water  and  the  range  of  X 
was  from  zero  up  to  about  0.95  and  n  ranged  from  0.005  up 

to  a 1 mo  st  1.0. 


A  -  3  Velocity  Pist ribut ion 

a.  Smooth  Boundary: 

The  velocity  data  on  the  centerline  of  the  channel  as  well 
as  off  the  centerline  were  plotted  with  u/u*  versus 
logy  U*.  ,  where  u*  is  the  shear  velocity,  equal  to  /r0/p  - 

at  the  foot  of  the  corresponding  normal  line,  and  P  is  the 
mass  density  of  water.  It  was  found  that  the  velocity 
distribution  agrees  well  with  the  well  known  Karman  - 

P  r and  1 1  equation: 


. 

- 
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=  5.75  1 o  g (  *  )  +  5*5 

u*  V 

even  for  values  of  X  as  large  as  0.95.  It  was  noticed  that 
near  the  f r ee-s ur f ac e  ,  there  is  a  dip  in  the  velocity 
prof  i  1 e  . 

b.  Rough  boundary: 

The  present  s emi-emp i r i c al  methods  for  computation  of 
turbulent  flow  in  open  channels  are  based  essentially  on  the 
work  of  Keulegan  (1938).  For  rough  boundaries,  Keulegan  uses 
the  Karman  -  Prandtl  velocity  distribution  equation  as 

u  y 

-  =  A  log(— )  +  C 

»*  k 

where  k  is  the  roughness  height  and  C  is  a  coefficient.  For 
the  densely  packed  sand  roughness  of  Nikuradse,  k  becomes  ks 
and  c-8.5.  Rajaratnam  used  an  averaged  value  of  u*  based  on 
y  r  S  where  Y  is  the  specific  weight  of  the  fluid,  r  is 
the  hydraulic  radius  and  So  is  the  bed  slope,  to  evaluate 
the  value  of  A.  He  found  for  4  types  of  roughness  the  value 
of  A  ranging  from  7.40  to  6.32. 


A_  4  ni  stribution  of  Boundary  She  ax  Stresses, 

Based  on  existing  data  from  literature  and  several 
supplementary  experiments,  a  plot  of  the  aspect  ratio  (B/D) 
versus  Tom/  Y  D  8„,  for  different  values  of  D/kg  in  shown  in 
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Figure  A-1  Centerline  Bed  Shear  Stress  for  Rectangular  Channels:  Smooth  &  Rough  Boundary 


Figure  A-2  Lateral  Distribution  of  Bed 
Shear  Stress  for  Rectangular  Channels: 
Smooth  Boundary 


Figure  A-3  Lateral  Distribution  of  Bed 
Shear  Stress  for  Rectangular  Channels: 
Rough  Boundary 
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Figure  A-l,  where  TQm  is  the  centerline  shear  stress.  From 

Figure  A-l,  it  could  he  seen  that  for  D/kg=3.0,  Tom/ Y  D  Sq 

becomes  equal  to  one  for  an  aspect  ratio  of  two  whereas  for 

the  channel  with  smooth  boundaries,  the  corresponding  value 

of  the  aspect  ratio  is  about  14.5.  In  the  case  of  the  smooth 

boundary  for  which  D/k  =co  ,  supplementary  experiments 

s 

confirmed  that  the  average  curve  obtained  by  Bajaratnam  and 

other  investigators  is  fairly  good.  Figure  A-2  shows  the 

distribution  of  T0/T0m  versus  on  the  bed  for  the  smooth 

boundary.  Figure  A-3  shows  the  distribution  of  T0/Tom  versus 

for  the  rough  bed  with  smooth  walls.  There  is  quite  a 

scatter  of  the  points  in  this  plot,  caused  basically  by  the 

variation  of  D/k  .  At  this  stage  of  the  work,  it  was  not 

s 

possible  to  label  the  different  curves  for  various  values  of 

D/k  .  An  average  curve  is  all  that  could  be  drawn, 
s 


£ _  5  Average  Velocity  Pred i ct ion 

For  smooth  boundaries,  u-(u*/k)  ln(y/yQ)  where  k=0.40, 
y0=m\)/u*  and  m=9.  If  U-average  velocity  in  the  vertical 
profile,  after  integration  with  proper  boundary  condition: 


u 


u 


* 


D  u* 

5.75  log (  — — 


+  3.0 


One  more  integration 
average  velocity  for 


with  U  from  z=0.0  to  z-B/2  will  give  an 
the  cross-section.  By  noticing  that: 


X  0  at  any  point  u* 

T  at  centerline  U*Q 

o  m 


and  fitting  an  equation 


to  the  curve  in  Figure  A-2,  one  can 
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proceed  with  the  integration  to  obtain  the  final  result  as 


V 


5.75 


log  ( 


D  u* 

v 


+  2.75 


( A  -  2  ) 


Compared  with  an  earlier  empirical  equation  (Chow,  1959) 
which  reads: 


V 


=  5.75  1 o  g ( 


r  v 


V 


)  +  2.75 


( A  -  3  ) 


wherein  v *  is  shear  velocity,  based  on  the  hydraulic  radius 
rather  than  the  depth  of  flow.  It  is  seen  that  the  constant 
3.25  is  higher  because  the  hydraulic  radius  r  is  always 

smaller  than  D ,  the  depth  of  flow. 

For  a  rough  boundary  the  governing  equation  is: 


-2—  =  5.75  log  (-2-)  +  C 

u*  k 

Upon  integration  from  the  bed  to  water  surface 
velocity  is  given  by: 


( A  -  4 ) 


the  average 


_U 

u 


-  =  5.75  log(* - )  +  6.0 

k 

*  s 


( A  -  5  ) 


(where  k  =Nikuradse  sand  roughness) 

once  more,  integration  from  the  centerline  to  the  wall, 
using  Figure  A-3  and  proceeding  in  the  same  way  as  in  the 
case  of  the  smooth  boundary,  the  average  velocity  for  the 
cross-section  is  given  by: 
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—  =5.75  log (-^-) 

u  .  k 

*o  s 


+  5.25 


(A-6) 


By  knowing  the  aspect  ratio  and  D/kg,  from  Figure  A  1  one 
can  obtain  TQm,  and  hence  u*q  and  using  Equation  A-6,  V  can 

be  calculated. 

These  results  are  similar  to  those  Keulegan  (1938)  who 
evaluated  these  coefficients  using  arbitrary  shear  stress 
profiles  and  the  experimental  data  of  Bazin. 


APPENDIX  B 


B-l  Channels  with  Flood  Plains  at  Two.  Levels. 


Often  successive  increments  of  overbank  deposition  are 
responsible  for  building  a  flood  plain.  The  erosion  of  the 
main  channel  bed  and  the  edge  of  the  flood  plain  creates  a 
narrow  strip  of  intermediate  flood  plain  level.  A  detailed 
study  by  Folman  &  Leopold  (1957)  showed  that  the  deposition 
on  the  flood  plain  does  not  continue  indefinitely.  The  flood 
plain  area  can  be  transformed  into  a  terraced  surface  by 
some  tectonic  or  climatic  changes.  A  terrace  overflows  less 
frequently  than  a  flood  plain.  When  a  flood  plain  and  its 
adjacent  terrace  are  flooded  during  a  very  high  flood,  the 
compound  cross-section  will  have  at  least  three  or  more 
different  bed  elevations* 

Only  one  experiment  was  conducted  with  the  multi-level 
flood  plain  channel.  The  dimensions  of  the  multi-flood  plain 
were  presented  in  Chapter  II.  (see  Figure  2. 2d).  Figure  B-l 
shows  the  distribution  of  u(y),  u(y')  and  u(y").  It  can  be 

seen  from  the  profiles  that  the  velocity  in  the  main  channel 
increased  with  y  and  the  maximum  velocity  will  be  obtained 
just  below  h  (first  flood  plain  elevation).  The  velocity 
will  continue  to  increase  at  lower  rates  up  to  y  -  h  and 
from  here  to  the  water  surface  will  remain  almost  constant. 
The  lateral  velocity  profiles  are  shown  in  Figure  B-2. 
considering  the  uppermost  flood  plain  (second  level  of  the 
flood  plain)  the  velocities  will  increase  continuously  as 
one  moves  towards  the  first  level  of  the  flood  plain.  The 
first  stage  flood  plain  evidently  acts  as  a  main  channel  for 
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Velocity  Profile  in  Main  Channel  ^ 


Figure  B-1 


Undisturbed  Velocity  Profiles:  Main  Channel 
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Figure  B-1  (cont.)  Disturbed  Velocity  Profiles:  Main  Channel 
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Figure  B-2  Lateral  Velocity  &  Bed  Shear  Stress  Distribution 
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the  second  stage  flood  plain.  The  velocity  continues  to 
increase  from  the  first  stage  flood  plain  to  the  mam 

ch  annel . 

The  variation  of  the  bed  shear  stress  with  2  ,  z'  and 

z"  is  also  shown  in  Figure  B-2.  In  the  main  channel  the  bed 
shear  stress  decreases  from  a  maximum  value  of  Tom  to  zero 
at  the  main  channel  wall  which  is  attached  to  the  first 
flood  plain.  In  the  first  flood  plain  channel,  Tc  decreases 
continuously  from  a  maximum,  at  the  junction  plane  to 
approach  an  undisturbed  value  of  Tq^j*  Then  T0  wil 
decrease  as  it  reaches  the  junction  plane  of  the  two  flood 
plains.  The  part  of  the  first  flood  plain  area  close  to  the 
main  channel  acts  as  a  flood  plain  for  the  main  channel.  The 
other  part  of  the  first  flood  plain  acts  as  main  channel  for 
the  second  flood  plain  level.  To  on  the  second  flood  plain 
level  will  continue  to  decrease  from  the  junction  plane  to 
approach  an  undisturbed  value  of  T0a)2  and  then  TQ  will 
decrease  as  the  bank  of  the  second  flood  plain  is 
a  p  p  r  o  a  ch  e  d  • 


p-2  Analysis  of  Experimental  Pesult?- 

From  lateral  velocity  distributions  one  can  see  that 

the  curve  consists  of  two  parts.  The  first  curve  looks  the 
same  as  the  lateral  velocity  distribution  of  a  main  channel 
with  a  flood  plain  channel  having  a  width  of  the  mam 
channel  equal  to  b'  and  the  depth  of  flow  equal  to  d'.  For 
the  second  part  of  the  curve,  one  should  consider  the  main^ 
channel  width  equal  to  b'  and  the  depth  of  flow  equal  to  d' 


,, 
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and  the  flood  flow  width  and  depth  the  same  as  the  second 
flood  plain  channel. 

The  lateral  velocity  distribution  curves  at  a  level 
y>(h+h')  for  the  above  two  single  compound  channels  could 
easily  be  calculated  and  plotted.  If  these  two  plots  are 
superimposed  in  such  a  way  that  the  origin  of  the  z  axis 
coincides  with  the  edge  of  the  first  and  second  level  flood 
plain  channels,  the  lateral  velocity  profile  of  the 
multi-level  flood  plains  will  form.  Using  the  results 
obtained  from  Chapter  IV,  the  two  separate  lateral 


velocities  were  calculated  and  plotted.  Superposition  of  the 
two  produced  the  final  lateral  velocity  distribution  which 
agreed  very  well  with  the  measurements. 


F -  3  Conclusions 

For  a  straight  multi-level  flood  plain  channel,  the 
interaction  between  any  pair  of  adjacent  channel  flows  will 
take  place.  The  channel  adjacent  to  the  main  channel  should 
be  treated  as  a  simple  flood  plain  for  the  main  channel  and 
the  same  channel  will  serve  as  a  main  channel  for  the  second 
flood  plain  and  so  on.  It  is  possible  to  obtain  the  lateral 


distribution  of  u  or  XQ  by  superposition  of  the  related 
curves  calculated  for  a  simple  compound  cross-section. 
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